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A SYSTEM AND METHOD FOR THE 
ANALYSIS AND PREDICT IO N OF ECONOMIC ^MARKETS 

FIELD OF THE INVENTION 

The present invention relates generally to a method 
and system for the analysis and prediction of economic 
markets. More particularly, the present invention forms a 
search space of different economic models, forms a fitness 
landscape by extracting observables from executions of the 
economic models with dimension reduction and by comparing the 
observables with actual data, and performs model selection 
and composite model synthesis through optimization over the 
fitness landscape. 



BACKGROUND 

The stock market and non financial markets are 
complex systems in which one or more stocks, goods, services, 
or their derivatives are sold or bought. Many techniques 
exist that attempt to analyze expected behavior of these 
markets, including predictions of price and volatility for a 
forward time series from any given date. Among these methods 
are those concerning prediction of volatility based on Black 
Scholes, and implied volatility given prices of options and 
the assumptions of Gaussian random walks in prices that are 
used to back estimate the volatility implied by actual option 
prices. In addition, procedures exist which posit some 
mathematical form for a feature such as volatility, such as 
truncated Levy, or Levy, or other distributions, then 
attempts to fit observed volatility distributions using these 
assumed underlying distributions and appropriate weighting of 
coefficients in the given mathematical form. 

Still other attempts to predict expected volatility 
uses Black Scholes formalism but treats variance in price as 
itself a random variable which changes over time, then 
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attempts to estimate the rate of change of this random 
variable. Still other efforts to predict volatility are 
based on causal models of various kinds. For example, some 
workers postulate a "herding" effect in which clusters of 
5 traders of various sizes form and trade in much the same way, 
causing small and large fluctuations in price, hence in 
volatility. Other attempts are based on Langevain equations 
which represent displacement of buy and sell orders at a 
moment from "equilibrium" , then seek stochastic equations 

10 describing the change of this displacement as a function of 
drift terms, price momentum terms, penalties for price 
variance, and volatility growth terms. Still other efforts 
are based on local linear models which embed a single stock's 
price series, or sometimes several price series, in a 

15 multidimensional space of the past N periods and attempts to 
locally fit the next price move of one or several prices. 

The predictive capability of each of these economic 
model can vary with the economic environment. Specifically, 
an economic model may exhibit a good predictive capability 

20 f° r ° ne set °f economic circumstances and a poor predictive 
capability for another set of economic circumstances. Since 
the economic environment changes with time, the utility of 
each economic model also changes with time. 

A great many researchers in economics, finance, and 

25 related areas have investigated different aspects of 

financial markets or markets in general. Topics of interest 
to the researchers and of relevance to our project have 
included diffusion of information in markets (Kyle, Albert S. 
(1985); "Continuous Auctions and Insider Trading," 

30 Econometrica, v. 57, 1315-35, Grossman, Sanford J., and 
Joseph E. Stiglitz (1980): "On the Impossibility of 
Inf ormationally Efficient Markets," The American Economic 

- 2 - 
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Review, v. 70,393-408), efficiency and statistical modeling 
of markets (Foley (1982) ) , speculation and bubble formation 
(Tirole, Jean (1982) : "On the Possibility of Speculation 
Under Rational Expectations, " Econometrica , v. 50, 1163-81. 
5 ) , and many others . 

Other existing work in market microstructure theory 
(O'Hara Blackwell, Maureen, Market: Microstructure Theory, 
Cambridge, MA, 1995 .) suffers from a lack of means for 
experimental verification (with some important exceptions, 

10 such as (Rust, John, John Miner, and Richard Palmer (1993) : 
"Behavior of Trading Automata in a computerized Double 
Auction Market", in The Double Auction Market) . For example, 
Glosten-Milgrom (GM) model examines the relationship between 
Dealer spreads and the proportion of informed Investors in 

15 the market- But it is very difficult to quantify the 

fundamental concept of "proportion of informed investors in 
the market" Ideally, one would prefer a numerical and 
experimental framework in which one could examine the regions 
of validity of analytic models like the GM . 

20 Another important shortcoming of the existing 

literature is the static nature of most of the models, i.e., 
results even for dynamic models are usually derived using 
some type of equilibrium assumptions. Such assumptions 
greatly limit one's ability to investigate dynamic, 

25 potentially non-equilibrium behavior of the system, and also 
one's ability to model realistically the market participants 
and as the result the market itself. 

Accordingly, there exists a need for a broad 
encompassing algorithmic framework to develop effective and 

3Q routinely improvable models and predictions of economic 

markets including stock markets and non-securities markets. 
In addition, there exists a need for a broad and improvable 

- 3 - 



SUBSTITUTE SHEET (RULE 26) 



WO 00/38094 



PCT/US99/30640 



algorithmic framework to seek means to price securities such 
as options and other derivatives. 

There exists a further need for a system for 
performing analysis and making predictions in a securities 
market using an agent -based models of the dealers and 
investors . 



- 4 - 
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SUMMARY OF THE INVENTION 



The present invention presents a system and method 
of economic analysis and prediction which dynamically adapts 
5 to a changing economic environment by selecting or 

synthesizing an economic model from a set of economic models 
based on the selected model's ability to make accurate 
predictions about an actual economic market. Specifically, 
the method and system of the present invention forms a space 
10 of different economic models, forms a fitness landscape by 
extracting observables from executions of the economic 
models, and performs model selection and composite model 
synthesis through optimization over the fitness landscape. 
The present invention includes a method for 
15 performing analysis and prediction of economic markets 
comprising the steps of: 

initializing a space of economic models ; 
executing one or more of said economic models to 
generate a corresponding one or more predicted time series; 
2Q determining one or more fitnesses corresponding to 

said one or more economic models by comparing said plurality 
of predicted time series to actual data; and 

searching for at least one of said economic models 

having a maximal fitness. 

The present invention further includes a system for 
performing analysis and making predictions of a securities 

market comprising 

a model of the securities market comprising: 

a plurality of dealer agents representing a 
30 corresponding plurality of market makers in the securities 
market, each of said dealer agents having at least one of a 



25 
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plurality of dealer strategies for setting a price of 
securities ; 

a plurality of investor agents representing a 
corresponding plurality of investors in the securities 
market, wherein said investor agent determines whether to buy 
or sell one or more of the securities at a given price; 

at least one market model comprising one or 
more rules for exchanging the securities; and 

a simulator for executing said model to perform 
analysis on the securities market . 

The present invention further includes a method for 
creating a portfolio of assets comprising the steps of: 

initializing ac lease one current portfolio of 
assets having one or more desired properties ; 

generating at least one next portfolio of assets to 
form said current portfolio, said next portfolio having a 
smaller size then said current portfolio; 

repeating said generating step after replacing said 
current portfolio with said next portfolio while said next 
0 portfolio has the one or more desired properties ; 

defining a size range having an upper limit equal 
to the size of said current portfolio and a lower limit equal 
to the size of said next portfolio; and 

searching for an optimal portfolio of assets to 

5 within said size range. 

The present invention further includes a method for 
optimizing a portfolio of assets comprising the steps of: 

determining a fitness landscape representation over 

a space of portfolios; 
Q determining at least one optimal searching distance 

in said fitness landscape representation; and 
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searching for optional ones of said portfolio of 
assets at said optimal distances. 

The present invention further includes a method for 
creating a portfolio of a plurality of assets comprising the 
5 steps of : 

determining a plurality of ant i - correlated families 
wherein each of said anti-correlated families contains two or 
more of said plurality of assets that are anti-correlated; 
and 

10 determining one or a more perspective portfolios of 

said plurality of assets from said plurality of anti- 
correlated families. 

The present invention further includes a method of 
modeling a portfolio of assets comprising the step of: 

15 defining a probability distribution of the value of 

the portfolio of assets with two or more Gaussian 
distributions. 



2Q BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 provides a data flow diagram describing a 
framework of the system and method of economic analysis and 
prediction of the present invention. 

FIG 2 shows a screen shot of the simulation of the 

2. b 

dealer-mediated securities market model. 

FIG. 3 shows the behavior of the population of 
Basic Dealers for different tick sizes. 

FIG. 4 displays the results of the dealer mediated 
30 securities market model for a market consisting of 7 Basic 
Dealers and one parasite. 
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FIG. 5 displays the results of. the dealer mediated 
securities market model for a market with a mixed population 
of dealers. 

FIG. 6 shows graphs of RMS (bid - True Value) and 
5 RMS (ask - True Value) as a function of the tick size for 

different populations of Market Makers in the dealer-mediated 
securities market model. 

FIG. 7 shows the average value of the spread as a 
function of the tick size for different combinations of the 
10 market makers in the dealer-mediated securities market model. 

FIG. 8 shows a graph of the standard deviation of 
the spread size as a function of the tick size for different 
populations of Market Makers m the dealer-mediated 
securities market model . 
15 FIG. 9 displays the main simulation control window 

of the dealer-mediated securities market model. 

FIG. 10 displays the Options menu of the dealer- 
mediated securities market model. 

FIG. 11 displays the Set Human Dealer Prices panel 
20 and the Set Tick Size panel of the dealer-mediated securities 

market model . 

FIG. 12 displays the "Show True Value" panel and 
the "trade history" panel of the dealer-mediated securities 

market model . 

FIG. 13 displays the "graphs menu" of the dealer- 

2 5 

mediated securities market model. 

FIG. 14 displays the "cash bar graph" panel of the 
dealer mediated securities market model. 

FIG. 15 displays the "Inventory Bar Graph" panel 
30 of the dealer-mediated securities market model. 

FIG. 16 displays the "Supply and Demand Graph" 
panel of the dealer-mediated securities market model. 
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FIG. 17 displays the U A11 Dealer Cash" graph of the 
dealer-mediated securities market model. 

FIG. 18 displays the "Dealer Assets" panel of the 
dealer-mediated securities market model. 
5 FIG. 19 shows che Dealer Inventory graph of the 

dealer-mediated securities market model. 

FIG . 2 0 displays a panel showing how the dealers 
track the true value reasonably well in the dealer-mediated 
securities market model . 
10 FIG. 21 displays a panel showing the simulation 

results for a smaller tick size in the dealer-mediated 
securities market model. 

FIG. 22 displays a panel of results obtained by 
executing the dealer mediated securities market model. 
15 FIG. 23 provides a flow diagram of a method for 

evaluating and minimizing risk. 

FIG. 24 provides the results of the method for 
evaluating and minimizing risk from executing on 500,000 
random draws from this "toy" world. 
2Q fig. 25 displays the histograms that determine the 

number of children to draw from during execution of the 
method for evaluating and minimizing risk. 

FIG. 26 displays a flowchart illustrating the 
method for portfolio modification. 
25 FIGs. 27, 28a and 28b display families of anti- 

correlated stocks that were generated by a method to create a 
portfolio of a plurality of assets with minimal risk. 

FIG. 29 displays a picture illustrating the anti- 
correlation relationship among family members that was 
30 generated by a method to create a portfolio of a plurality of 
assets with minimal risk. 
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FIG . 3 0 discloses a representative computer system 
3010 in conjunction with which the embodiments of the present 
invention may be implemented. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

System Framework 

FIG. 1 provides a data flow diagram describing an 
overview of the system and method of economic analysis and 
prediction of the present invention. In step 110, the method 
executes each economic model 105 in a set of economic models 
to form a corresponding family of predicted time series of 
desired features 115. Without limitation, the economic 
models 105 include a variety of different model types such as 
parameterized models, agent-based causal models, aggregate 
causal models, and intermediate causal models. The economic 
models 105 further include models based on a variety of data 
mining tools such as time series analysis. The economic 
models 105 also include causal models based on features such 
as the economic web model in which N stocks are located, 
hence the sensitivity of the values of the stock to the 
introduction of new goods and services and the extinction of 
old ones as described in copending patent application, 
-System and Method for the Synthesis of an Economic Web and 
the Identification of New Market Niches", U.S. Application 
Number 09/080,040, filed 5/15/98, the contents of which are 
herein incorporated by reference. The economic web models 
are further described in "A System and Method for 
Coordinating Economic Activities Within and Between Economic 
, Agents", U.S. Application Number 60/091,656 and International 
Application number PCT/US99/15096 , the contents of which are 
herein incorporated by reference. 

- 10 - 
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The economic models 105 further include causal 
models based on detailed agent based models of the "players" 
in a stock market, such as dealers, traders, plus their order 
books, plans and facts concerning wishes to buy or sell 
5 stocks or derivatives as characterized by a finite list of 
features including price, volume, strike price, dates, 
underlying value of the stock by various estimates, and other 
factors. In addition, these agent based models may include 
strategy spaces for the different agents, means by which 
10 agents change strategies, and may include sets of hypotheses 
that agents hold with respect to other agents prospective 
actions, other agent's strategies, and other agents 
hypotheses. . Specific examples of economic models 105 
include AR(n), ARMA, ARV , ARCH, GARCH, BlackScholes , 
15 Langevain equations, etc. Accordingly, the space of economic 
models 105 is broad and indefinitely extendable. 

The execution of step 110 yields time series data 
for models 105 having the same types and different parameter 
values, models having different types, etc. For the agent 
20 based models in the set of economic models 105, step 110 

executes on models with different strategy spaces, different 
inventory requirements, different preferences, etc. For the 
agent based models where the agents have different models of 
each other, step 110 executes with different sets of rules 
25 which define how the models of the agents change over time. 
Since the models change with time, step 110 executes the set 
of economic models 105 for different time periods. For 
example, step 110 could be executed for a time period 
beginning 10 years ago and running until the present time. 
30 Likewise, step 110 could also be executed for a time period 
beginning 1 year ago and running to the present time. 
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In step 120, the method extracts a set of 
observables from the predicted time series and corresponding 
time series of the actual economic market to create a fitness 
landscape 125 of the space of economic models 105. Execution 
5 of step 120 yields metric distances between different models 
for different time periods. Specifically, execution of step 
120 yields metric distances between models having the same 
type and different parameter values, models with different 
strategy spaces, models with different inventory 

10 requirements, different preferences, etc. 

The method extracts a set of observables from the 
families of predicted time series and creates a fitness 
landscape using a wide variety of data mining techniques and 
dimension reduction techniques. Dimension reduction 

15 techniques include principal component analysis, localized 
principal component analysis, independent component analysis 
and the sparse basis of landscapes as described in copending 
patent application, "A System and Method for Coordinating 
Economic Activities Within and Between Economic Agents", U.S. 

20 Application Number 60/091,656 and International Application 
Number PCT/US99/15096 , the contents of which are herein 
incorporated by reference. 

In an alternate embodiment, step 120 computes 
probability distributions for the observables for the 

25 families of predicted time series. In this embodiment, step 
120 computes a comparison between different probability 
distributions to approximate metric distances between 
corresponding predicted time series. Preferably, the Randan- 
Nikodyn Thereom is used to approximate the metric distance. 

3Q Alternatively, the Kulback-Liebler method could be used. 

However, as is understood by persons of ordinary skill in the 
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art, other techniques could be used to compare probability 
distributions and approximate their metric distances. 

Execution of step 120 assesses a given state of an 
economic market. Without limitation, that current state 
might be an embedding of the current and N past day's prices 
of M stocks in an (N+1)M dimensional state space showing by a 
single point the entire history of all M prices over the past 
N +1 periods. The current state of the market, however, may 
be broader, and include prices, volumes, and other desired 
features of the stocks or goods sold and bought. The current 
state might be the present state, or a past historical state 
or a proposed future hypothetical or expected state of the 
market. More generally, the embedding might include the N+l 
past moments and F future moments or periods, as sampled from 
historical data such that the state space has (N+l+F) * M 
dimensions . 

Again at the highest level, the chosen model is 
utilized, together with an assessment of the current state of 
the market, to predict one or more successive states of the 

0 market. This process may be used a single time to predict, 
for example, a time series of the M stocks for P succeeding 
periods, where P is 1 or more, and periods are chosen as 
desired for past and future analysis. 

More generally, the predicted future time series 

5 might be reformed deterministically , or stochastically, from 
the current state based on features of the model's mapping of 
past to future states. Thus, for a local linear method, the 
statistical features of a mapping of a set of neighboring 
past points to a set of future points in the (N + 1 + F) M 

0 state space is used to pick, stochastically predicted future 
points from the current state to, generally, a plurality of 
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different members of "walks" from the current state into the 
future . 

The ensemble of these walks may then be used to 
measure features such as price. More particularly, the 
5 ensemble of these walks, or partial differential equations 
based on the cumulants of the ensemble can be used to 
characterize the ensemble. Based on numerical study, or the 
partial differential equations, the ensemble can be used to 
assess hedging, risk minimization to a single option, or a 
10 portfolio of holdings. 

In addition to the use of one or an ensemble of 
predicted future price series, or time series of other 
features such as volume and otherwise to assess risk and 
proper prices for options and other instruments, the 
15 predicted time series can be compared to the observed 

behavior of the real market 118 in order to test, modify and 
improve the initial choice of model for the market. 

At the highest level, this testing procedure 
consists in the use of one or more comparisons between the 
20 predicted time series 115 and actual time series 118, as 
assessed by features as simple as direct visual comparison, 
to principal component analysis, independent component 
analysis and other sophisticated modes of comparison of data 
sets. These comparisons yield one or more "errors", that 
25 measure different aspects of the deviation of the predicted 
behavior of the market to the observed behavior of the 
market . 

Given a set of errors, the highest level algorithm 
constructs a "fitness measure" of the model in step 120. This 
30 fitness measure may, in general, consist in a single number, 
or might be a vector with all of the different features of 
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errors noted by the comparison of fact and model predictions, 
or a still more general measure of fitness of the model. 

Again at this highest level, the fitness measure or 
measures creates a fitness landscape 125 over the set of 
5 models 105 of the market. Given any one of a plurality of 

means of defining "nearby models" , the most general algorithm 
achieves a space of models of the market. Since each model 
has a fitness, given the space of models 105, the algorithm 
achieves a fitness landscape 125 over the space of models 

10 105. in this fitness landscape 125, it is natural to 

represent minimization of errors by being fitter, and allow 
the error measures to be a height . Hence, in this 
formulation, "fitter" is lower, corresponding to lower 
errors. Without loss of generality, fitter can be taken as 

15 higher by reformulating the error measures. Without loss of 
generality, let "high" be fitter than "low", so peaks 
represent high fitness and thus, low error. Then the fitness 
landscape 125 may, in general, be smooth and single peaked, 
rugged and many peaked but correlated, with correlations 

20 lengths that differ in different directions, or random/ 

Means to search a landscape, to characterize and to represent 
a landscape are described in copending patent application, 
"A System and Method for Coordinating Economic Activities 
within and Between Economic Agents", U.S. Application Number 

25 60/091,656 and International Application Number 

PCT/US99/15096 , the contents of which are herein incorporated 
by reference . 

In step 130, at the highest level, the algorithm 
utilizes, first, means to assess the statistical structure of 

30 the fitness landscape 125 over the space of models 105 of the 
market, such as its correlation lengths in different 
directions. The algorithm determines the optimal search 
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distance as a function of the current fitness of a model of 
the market and the correlation structure of the landscape, 
etc. as described in "A Method for Optimal Search on a 
Technology landscape", International Application 
5 PCT/US99/19916, filed August 31, 1999, the contents of which 
are herein incorporated by reference. Then, based on these 
features, the highest level algorithm chooses, by mutation, 
recombination, or other methods, a new trial model of the 
market to test and assess the fitness of that model. 
10 in step 130, over a succession of trials, the 

highest level algorithm searches for fitter models of the 
market. In general, this search procedure may be a local 
walk via fitter neighbor variants that becomes trapped on a 
local fitness peak in the space of models of the market or a 
15 more general search procedure that persistently searches in a 
space of models, generates new models, and continues to 
search for good models. This persistent search in the space 
of models in parr reflects the fact that the market itself 
may be persistently non_stat ionary . For example, its rules 
2 0 and behaviors may be changing. Thus, the highest level 

algorithm can tune the rate of search in the space of models 
to match the non_stat lonarity of the market itself such that 
the search procedure persistently finds good models of the 
market . 

25 in addition, at this highest level, the economic 

models 105 may be algorithmically augmented by a variety of 
algorithmic procedures such as changes in the coefficients of 
a model or changes in the functional mathematical form of a 
model. For example, step 13 0 can alter a term in an equation 

30 from linear to squared, or generate entirely new models by 
mutation and recombination of terms and elements or agents 
that are present in existing models. 
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In step 130, the method searches the space of 
economic models and synthesizes new economic models 135. The 
search and synthesis step 130 will also include the use of a 
stacking function. The stacking functions can be arbitrary 
5 functions of the models of the system. Genetic operators such 
as recombination and mutation operate over the space of 
stacking functions to produce new composite economic models 
135. The stacking functions must have adequate diversity to 
span the dimensionality of the behavior space. 

10 The scope of some of the searches from step 130 

will be restricted to intra-class searches. These searches 
include the refinement of parameters for parameter type 
models by comparing the actual data 118 with the families of 
predicted time series data generated by the model. In other 

15 words, the fitness of each economic model 105 is the metric 
distance of its predicted time series 115 from the actual 

time series 118. 

The search and synthesis step 130 will also include 
the use of biological functions. The use of biological 

20 functions produces hybrid economic models 105 which are 
composed of individual economic models. Each individual 
component economic models 105 contributes a portion of the 
predictive time series data produced by the hybrid economic 
model 105. The biological functions adjusts the permutation 

25 and percentage contribution of each component economic model 
105 in the hybrid economic model 105. 

The extraction of observables from ensembles of 
time series generated from hybrid economic models 105 and 
composite economic models 105 can be done in different 

30 manners. In one exemplary technique, the method takes a 

weighted sum of each of the component economic models 105 and 
subsequently performs extraction of the observables. In an 



- 17 - 



SUBSTITUTE SHEET (RULE 26) 



WO 00/38094 



PCT/US99/30640 



alternate technique, the method performs extraction of the 
observables for each component economic model 105 and 
subsequently takes of weighted sum of these sets of 
observables . 

5 The non-stationarity of economic markets will 

induce a flow in the behavior space created by step 120. 
Specifically, the values of the observables extracted through 
the data mining and dimension reduction process will change 
with the economic markets. The present invention can be used 
10 to analyze the non-stationarity of the economic markets by 
studying the flow in the behavioral space. In turn, the 
present invention can be used to predict the non-stationarity 
of the economic markets from the observed flow in the 
behavioral space. 

15 

Algorithmic Procedures for Options, Hedging and Other 
Investment Decisions 

Different economic models 105 can be used with the 
framework of the system and method of the invention shown in 

20 FIG. 1. A set of past time series and matching statistics 
forms an exemplary economic model 105, This process begins, 
without limitation, with a time series at some fixed 
interval, say a day, or an hour, or a week, with the prices 
of one or more assets, say of N different stocks. For each 

25 stock, there exists the past time values at the current price 

at time T, the price at T - I , T -2, .T -3 T - tau. In 

addition, there exists the future time prices T + 1, T + 2, T 
+ 3, . . . . T + tau. In addition to merely including price 
information at these different times, it is possible to 

30 include further information such as volume and price per 
trade at these different times, or lumped into small 
intervals of time "boxes". In addition, it is possible to 
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include other features such as measures . of expectations of 
changes in interest rate, or the deviation of an index such 
as the DJIA above or below its average value for some past 
period or periods. It is possible in addition to include 
measures of concern over political events, economic events 



etc 



Without limitation, and for ease of description 
only, we now consider the case where only the prices of N 
assets for the past and future 2 x tau + 1 time moments are 
considered. This data is obtained from the actual historical 
data record for different moments, T, chosen, for example on 
different days in the past year or several years or more. The 
accumulation of these data yield a mathematical state space 
structure which has a -past" a -now" and a -future". The 
dimensionality of this state space is N (2xtau + 1) . Thus, if 
tau is chosen to be 4 , meaning using data on prices from the 
moment T, back to T - 1 , T -2 , T -3 , T-4 and forward to T +1 , 
...T+ 4, then each asset corresponds to 9 prices, or 9 
variables, and the total state space has 9N dimensions. 
0 Accumulation of data over a past set of times,- say 

the past several years at daily intervals, will yield a set 
of points, one for each moment in the past chosen as -now". 
In general, nearby -past" points will have more or less 
nearby "future" points. Thus these data clouds of lumped 
5 nearby past and future points convey some information about 
the mapping from past to future in this state space. 

In particular, the mapping from a nearby set of 
-past" points to their corresponding -future" points will be 
used in the present algorithm to generate a stochastic 
0 process that, at each step from a -now" point, matches the 
local statistics of the past - future mapping. 
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In general, the number of sampled points in this 
state space will be sparse. Interpolation procedures will be 
useful for the present algorithm which will yield good 
approximations to the statistics of "past" - "future" 
5 mappings for points in the state space that have not been 
sampled. In general, markov random fields and other 
procedures introduced below allow this interpolation before 
or during the stochastic process next constructed from this" 
"past" , now, future high dimensional state space. 
10 The next step in the algorithm begins at any chosen 

date, in particular, the actual current historical date, and 
utilizes the state space to generate a "walk" beginning at 
the actual or chosen "now" , and choosing a "next price" 
randomly to fit the local statistics of "past", "Now, 
15 "future" prices in its vicinity of the state space. This 

random choice moves generates at next "mock" price at mock T 
+ 1, hence at a new point in the state space corresponding to 
the new price, and the new change in price, i.e. price at T 
price at T - 1, price at . . . . T - 2 , . . . T - tau. Thus, the new 
20 mock price generates not only a new price, but locates the 
stochastic process at a new point in the "past" of the state 
space . 

From the new point in the past of the state space, 
the stochastic process chooses at random, a mock future price 

25 matching the statistics of the future from this new point in 
the past xn the state space. This stochastic process is 
iterated as many times as desired to yield a mock time series 
of "future" prices of the set of N assets from the initial 
position, "now" in the state space. 

3Q The next step in the algorithm generates a second 

"mock time series" from the same initial position, "now" in 
the state space. The algorithm then generates as many mock 
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time series from the same initial position, "now" in the 

state space. 

The family of mock time series forms the basis for 
5 all further evaluations of investment and hedging strategies, 
as well as measures of volatility, and means to obtain 
options prices. For example, using this family of mock time 
series, coupled with an assumed option, date of sale and 
strike price, one can numerically investigate the expected 
10 return on that option for each member of the mock time 

series, hence can derive the average and the higher moments 
of the return on the option over the family of mock time 
series. Since there are a finite, but indefinitely large 
number, of members of the mock family of time series, it is 
15 always possible to obtain these expected values on the mean, 
variance, and higher moments. 

More generally, the motions in the state space 
exhibit the correlations in price and price changes among the 
N assets. Thus the family of mock time series can be used to 
20 construct good hedging strategies. More, the use of genetic 
algorithms and other procedures known in the art can be used 
to choose optimal hedging strategies among the N assets and 
their derivatives, given performance criteria such as risk 
aversion, and other desired features of the instrument. 
25 Further, in general, the real markets are non- 

stationary. Thus, in general, over real time the mapping 
known "Past" to "Future" of the state space changes. These 
changes can be thought of as a higher order map that 
transforms the "Future" of the state space onto itself. This 
30 higher order map can be created as a Tensor process mapping 
the "Future" state space onto itself. In turn, the Tensor 
mapping may not be stationary, and a still higher order 
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mapping of the Tensor onto itself can be created. These 
Tensor-Tensor maps are then used to modify the state space 
"over time" from the "initial now" to generate an improved 
family of mock/time series. 
5 Indeed, the first construction of the state spare 

used historical data over, for example, the past several 
years. The market mapping of "past to future" is, in 
general, non- stationary over this historical past. Nearby 
points in the "past" in the state space correspond, in 
10 general to real time points that may be separated by months 
or years. Thus the "dispersion" of "Future" data points may 
reflect randomness of the market at any one time, or the non- 
stationarity over time. 

The "Tensor" method of mapping the "future" of the 
15 state space onto itself can be used to analyze and improve 
the initial state space by partitioning the state space into 
real time segments, each having a modest length, (i.e. 1 
month compared to several years - say 10 years) . Then the 12 
x 10 - 120 "past-future" state space maps can be compared via 
20 the "Future - Future" tensor matching procedure to find; 

estimate and modify the initial "past - now - future" state 
space to achieve an improved "Tensor" mapping and higher 
order mapping to use to generate families of mock time 
series upon which asset strategies and pricing can be based. 

25 

Modeling and Simulation 



The space of economic models 105 of the present 
invention goes far beyond simply fitting time series to 
stochastic models. Such approaches, including AR [n] , ARV, 
ARMA, ARCH, GARCH , LSSM, Levy, and truncated Levy models, 
merely fit the data (usually of a stock price, or of price 
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differences) to a mathematical model. Then that model is 
used in various ways for volatility prediction, risk 
assessment, option pricing, and optimal hedging. But those 
models generally have little or no basis in how markets 

5 really work, and treat the price data merely as a 
mathematical time series to be analyzed and fitted. 

In contrast, causal models include ideas, 
re p resen ted in equations, about how markets really work, and 
about how agents in those markets make decisions . The 

0 observed data, such as the price as a function of time, is 
then understood in terms of its underlying causes, not just 
as a set of numbers. By refining and calibrating our causal 
models against real-world information and data on market and 
agent behavior, we can obtain a much better understanding of 

5 price series than models without any causal mechanisms. This 
better understanding in turn provides more predictive power, 
and thus improved volatility prediction, risk assessment, 
option pricing, and optimal hedging. 

Our causal models of markets span three levels of 

0 modeling detail: 

1. Individual Agent Level 

2 . Agent Distribution Level 

3 . Aggregate Level 

These are discussed in more detail below. The present 
5 invention uses models at all three levels. The present 

invention calibrates models at levels 1 and 2 against each 

other, and similarly at levels 2 and 3. 

The division into three levels is actually somewhat 

arbitrary. In particular, there are many possible levels of 
0 detail within level 2, and many hybrid cases between levels 2 

and 3 in which some variables are modeled as distributions 

and some as aggregates. 
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The individual agent level models represent (in the 
computer) individual agents who buy and sell securities. 
These agents can include investors, brokers, specialists, and 
market makers . The security prices are determined by the 
5 activity of the agents, either through specific specialist or 
market -maker strategies, or through tatonnment rules in which 
the price of a security moves upwards when demand exceeds 
supply and downwards when supply exceeds demand. 

The agents each have a current portfolio (or 

10 "inventory") of holdings in each of a number of financial 
securities (stocks, bonds, options ,- etc) , and also have 
available reserves of liquid assets ("cash"). They change 
their portfolio, buying and selling securities, according to 
a current "strategy". These strategies may be represented by 

15 parameterized rules, a neural network, or fuzzy logic. The 

strategies, however represented, may involve agents' internal 
predictions of the future behavior of the market and/or of 

other agents. 

Each agent may contain a set of many different 

2 0 strategies at the same time, using those that seem, on the 
basis of past experience, best under any particular 
circumstances. Alternatively a weighted average of the set of 
strategies may be used, weighting each according to its 
observed past success and generating composite or 

25 intermediate actions. 

In general the set of strategies, and the 
individual strategies within the set, change through time as 
internal models are refined through learning, and as 
different approaches are found to be more or less successful. 

30 Entirely new strategies are generated using genetic algorithm 
and similar evolutionary programming techniques. The new 
strategies replace old unsuccessful ones. 
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At the opposite extreme, the least detailed level, 
the models at the aggregate level compress the behavior of 
the agents into one or a few aggregate variables. There is 
no representation of individual agents, but their average or 
5 aggregate behavior is represented. These are much simpler 
mathematically than the individual agent models. These 
models are further described in R. Cont and J. -P. Bouchaud, 
Herd Behavior and Aggregate Fluctuations in Financial 
Markets, preprint cond-matl97123i8 , available from 
10 http://xxx.J-anl.gov, 1998, the contents of which are herein 
incorporated by reference. The models at the aggregate level 
are further described in J. -P. Bouchaud and R. Cont, A 
Langevin Approach to Stock Market Fluctuations and Crashes, 
preprint cond-mat/9801279 , available from 
15 http://xxx.lanl.gov, 1998, the contents of which are herein 
incorporated by reference. 

The models at the aggregate level usually consist 
of a small number of coupled nonlinear differential or 
difference equations, some of which involve stochastic 
20 forces. Overall, these are Langevin- type equations in the 
terminology of physics. Some terms xn these equations may 
also involve integrals or averages over past behavior, 
usually in the form of weighted moving averages. Analytical 
methods of the Ito calculus, including construction of 
25 Fokker-Planck equations, may be applied in some cases, but in 
general numerical solutions are used. 

The aggregate variables that appear in the equations may 
include (but are not limited to) : 

1. aggregate demand- the total number of shares that 
3Q the trader agents want to buy, minus the total 

number that they want to sell. 
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2. inventory level-the total amount of inventory 
(counting short sales as negative) that the market- 
maker or specialist agents currently hold. 

3. average investor confidence-a "psychological 11 
variable depending on the recent stability or 
predictability of the market, and affecting the 
agents' risk aversion (their balance between desire 
for expected profit and avoidance of loss) . 
Similar variables may also be used for market- 
makers or specialists. 

4> average net capital entry rate -the average rate at 
which new capital enters the market, less the exit 
rate. 

Besides these and other agent -aggregate variables, there are 
of course pure market variables, including price, price 
trend, price volatility, price kurtosis, price correlation 
time, and the corresponding quantities for price differences 
and squared prices. There are also adjustable parameters in 
the model, including parameters to represent average growth 
rates, average shock and noise rates, strengths of couplings 
between various pairs of variables, etc. 

Together the coupled equations produce a causal 
model of the agents' aggregate behavior and the market's 
response. In fact each drives the other; the market reflects 
the behavior of the agents, and the agents take their cues 
from the market . 

At the intermediate level of detail of the models 
at the distribution level, the fundamental variables are 
distributions of aggregate quantities, such as those listed 
above for the aggregate level. These distributions are used 
in two ways : 
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1. to represent the variation of a quantity (such as 
average investor confidence) or a parameter (such 
as risk aversion, or trading lot size) across 
agents ; 

5 2 . to capture the uncertainty of such a quantity or 

parameter through time, in the form of a 
probability distribution. 

By using distributions rather than single aggregate 
quantities, we can model real agents and markets much more 

10 realistically, without needing to go all the way to the 

individual agent level. It seems indeed almost essential to 
have a distribution of large and small trades in order to fit 
the observed non-Gaussian distribution of stock price 
changes. Having a fixed size leads (according to the central 

15 limit theorem) to Gaussian price changes unless "herd" 
effects are invoked as described in R. Cont and J. -P. 
Bouchaud, Herd Behavior and Aggregate Fluctuations in 
Financial Markets, preprint cond-mat 19712318 , available from 
http : //xxx. J. anl . gov, 1998 (the contents of which are herein 

20 incorporated by reference) , and a herd acting together is 
effectively creating a large trading block in any case. 

The equations for the distributions are derived 
from the corresponding Langevin- type equations for the 
variables whose distributions we are . representing . This 

25 leads to coupled differential or difference Master Equations, 
with convolution integrals involved in the coupling between 
variables. In most cases, the present invention replaces 
these Master Equations by coupled differential or difference 
equations for the cumulant moments of the distributions 

3 0 involved, truncating the cumulant expansion after the first 
few terms (keeping at least 4 terms, up to at least the 
kurtosis) . This procedure leads to algebraic evolution 
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equations for the cumulants, enabling us to predict, for 
example, the price distribution of a stock. 

At the individual agent level, and at the aggregate 
level, the causal models produce time series of securities 
5 prices as a result of agent behavior and trading. 

Preferably, the models are run many times in order to obtain 
an ensemble of such time series. This ensemble may be 
regarded as a sample from an underlying probability 
distribution for the whole stochastic process (the causal 
10 model) that generates it. The present invention usually fits 
such an ensemble of time series to a stochastic process 
model, such as (but not limited to) a truncated Levy flight. 
Neural network techniques, including (but not limited to) the 
expectation-maximization (EM) algorithm and Bayes networks 
15 may be used to perform the fitting. 

In contrast, the intermediate level approach, 
working directly with distributions, generates a probability 
distribution for the stochastic process directly. 
In all cases such a distribution is fundamentally what is 
20 required for accurate volatility prediction, risk assessment, 
option pricing, and optimal hedging. 

The time-series ensembles, or their stochastic- 
process fits, are compared with real data from actual 
markets. The statistical differences between the model and 
25 the real data are then used to refine the model parameters or 
change the architecture of the model itself. 

A Dealer-Mediated Securities Market Model 

The space of economic models 105 further includes a 
30 dealer-mediated securities market model. The advantage of 

the agent -based models of this model is that a wider range of 
phenomena will become open to both investigation and 
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validation. For instance, an agent -based model allows us to 
incorporate learning on the part of Dealers and Investors,, as 
well as to investigate dynamics, which otherwise may not be 
amenable to existing analytical techniques. Particularly 
5 when individual agents learn or adapt over time, the 

resulting co-evolution is normally completely intractable. 
' The computer simulation becomes the scientific tool by which 
we may model, and begin to understand, such systems. 

Preferably, this market model is implemented with a 

10 platform in the Java programming language for simulation and 
investigation of such a market. The simulation consists of a 
computational Agent -Based Model (ABM) of a market that trades 
securities. We have used the model to investigate questions 
involving dealer strategies, especially parasitic strategies, 

15 and their implications for real -world issues for stock market 
exchanges such as Nasdaq. The present invention investigates 
decimalization (the reduction of the minimum quote increment 
or tick-size to It) and the market's effectiveness at 
performing price discovery under different proportions of 

20 parasitic strategies. Our results indicate that a move* 

towards decimalization may have a significant negative impact 
on price discovery in the presence of parasitic strategies, 
and perhaps slightly impede it in a non-parasitic market. 
The model can also be used to investigate other issues. 

25 These models, and the need to validate them 

experimentally, arise out of very practical considerations 
about the micros tructure of security exchanges such as 
Nasdaq. In recent years, Nasdaq has changed the way it 
works, e.g., the new order handling rules and the recent 

3Q proposal to decimalize the system by reducing the tick size 
to a penny. Every one of these structural changes has the 
potential to produce sweeping effects in the structure of the 
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market - consequences not necessarily intended by their 
implementers . 

To this end, the present invention includes an 
agent-based model of a dealer-based market. This computer 
5 simulation represents the players in the market, Dealers and 
Investors, as "agents." Agents are computational objects 
that behave so as to maximize some fitness criterion. In the 
case of a market, the fitness is (at zeroth order) 
profitability, with additional criteria, such as risk 
10 aversion, appearing as higher-order corrections. This model 
provides an "economic laboratory" to explore the behavior of 
the market and the consequences of modifying its institutions 
and rules . 

The main questions we address concern the effects 

15 of tick-size, parasitism and their interaction with the 

market rules and structures. In addition, these models are 
quite capable of providing research tools for many other 
topics. The models are particularly adapted to provide rea- 
sonably robust answers to the focal questions of our 

2 0 investigation . 

Parasitism is broadly defined as a strategy which 
attempts to make a profit without contributing to the process 
of price discovery. We will also use parasitism with a 
slightly different, albeit intuitively clear meaning - as an 

25 ability to take advantage of other players' actions. As 

discussed later, in the latter sense some players may behave 
as parasites, even without "consciously" realizing or 
attempting to do it where "consciously" means that its rules 
may not explicitly imply any parasitic behavior. An example 

30 of an explicitly parasitic strategy might be a strategy in 

which the trader hangs back (behaves non-competitively) until 
market conditions are ripe for exploitation and then trades 
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inside the spread, making a quick profit. In the context of 
double-oral auction markets, such strategies were observed to 
be highly effective, outperforming all others as indicated in 
Rust, John, John Miner, and Richard Palmer (1993) : "Behavior 

5 of Trading Automata in a computerized Double Auction Market", 
The Double Auction Market: Institutions, Theories and 
Evidence, Reading, Massachussett s : Addison-Wesley , the 
contents of which are herein incorporated by reference. 

The main intention of the model of the present 

0 invention is to create a population of interacting agents 
(Investors and Market Makers/Dealers) that together produce 
price discovery in the marketplace. The Java simulation we 
have implemented is very powerful, and can incorporate 
additional features of the real world, such as a limit-order 

,5 board, trading of multiple shares, , or trading in a wide range 

of volumes etc. 

The market contains a plurality of securities whose 
value is exogenously specified: the underlying value of the 
security is assumed to fluctuate according to a random walk 
0 with a drift. The market consists of a finite population of 
players trading in that single security. We have two types 
of players: Dealers (Market Makers) and Investors. 

Each Investor has access to the current "true" 
value of the underlying security adjusted by an idiosyncratic 
5 error, which in turn is determined by the investor's degree 
of inf ormedness . The simulations were run with errors 
distributed by a Gaussian law. The simulations can also 
investigate, information time lags as well as errors. 

Dealers, on the other hand, do not have any prior 
1Q information about the "true" price. To a certain extent, 
competition drives Dealers to discover the price through 
market interactions (this is the process of price discovery 



- 31 - 



SUBSTITUTE SHEET (RULE 26) 



WO 00/38094 



PCT/US99/30640 



which markets mediate) . Dealers are obliged at every moment 
to post their current bid and ask prices, and their volumes, 
on the public board. Investors decide whether they consider a 
purchase or sale of the security beneficial to them (by 
5 comparing their knowledge of the true value with the current 
market quotes), and, if so, they execute a trade at the best 
available price . 

The timing of the simulation is determined by a 
scheduling agent, which assigns to every player (dealer or 
10 investor) a period of time in which the player is able to 

take his next action. An action for a dealer is to choose his 
bid and ask price, and for an investor is to buy or sell at 
the best currently available price. A dealer can also update 
his quotes after every transaction he performs. There is 
15 much flexibility in the model with respect to the relative 
frequency with which each dealer and investor acts. In most 
simulation runs, dealers acted at a fast rate, with investors 
having a larger range of rates, from very fast to quite slow. 
The scheduling process could be adjusted to model the order 
20 flow of different kinds of security. 

Each player is endowed with a fixed quantity of 
money and a fixed initial inventory. The total number of 
units of the security can be fixed in the simulation. If at 
any time a player's total assets fall below a negative 
25 constant (fixed at. -$1000 in most of our runs), that player 
is removed from the simulation and replaced by another who's 
strategy parameterization is randomly selected. The old 
player's inventory (positive or negative) is redistributed to 
the remaining players (who are charged/compensated 
3 q accordingly) . Since the new player has an initial endowment 
of cash, money flows into the system as a consequence of 
player bankruptcy. 
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The model market consists of a framework in the 
Java programming language. This framework includes the model 
market, which is populated with Dealer and Investor agents 
who trade in one security; an agent scheduler, and a 
5 graphical user interface that allows one to control the 

microstructure of the market and monitor the progress of the 
simulation from both macroscopic (statistical) and 
microscopic (agent) points of view. The core simulation code 
can also be run separately from the graphical user interface 

10 (in "batch mode"). 

At the highest level, the simulation consists of 
four types of software objects: Price, Dealers, Investors, 
and the Market itself. At the next level down, all of these 
objects have member objects that provide communication and 

15 bookkeeping capabilities. A detailed breakdown of the 

simulation and its component objects is given later. FIG. 2 
shows a screen shot of the simulation of the dealer-mediated 
securities market model. Each menu and chart shown is 
described in detail later . 

20 The models represent agents' fitness functions and 

strategies . We have two types of strategies : Dealer 
strategies and Investor strategies. Both types of agents are 
concerned with making money. However, Dealers, as the 
mediators of all trading, have much more latitude in their 

25 strategies. Dealers must set both bid and offer prices and 
volumes. Investors must decide only on the course of action 
at any given moment, namely, whether to buy or sell at a 
given price and how many shares. 

The underlying value follows a stochastic dynamic. 

3 q At each opportunity to trade, an Investor is given this value 
plus some amount of error corresponding to his degree of 
"informedness." If a Dealer's current prices are below 
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(above) the Investor's perception of value, then the Investor 
will buy (sell) one or two shares of the security. Using this 
strategy, we can easily set the proportion of informed 
Investors, as in the GM model, while retaining a much more 
5 realistic model of market dynamics. The present invention 
also comprises investors who wish to trade a range of 
volumes . 

The core of our current work focuses on Dealer 
strategies, where there is a very rich environment for re- 

10 search. Our main goal was to investigate the effects of 
parasitism, tick size, and the combination of both on the 
market. For that, we have developed two broad classes of 
dealer strategy: one that makes price discovery happen, and 
another that uses information created by the rest of the 

15 population (parasitic strategies) . As discussed later, there 
is not a clear distinction between the two. With this goal 
in mind, we developed and implemented several specific Dealer 
strategies: 

• A "Basic Dealer" only reacts to trades which come his 
own way. He maintains his quotes until he receives a 
certain number of trades on the buy (or sell) side, and 
then adjusts the quote on that side appropriately. It is 
a simple strategy which manages to track price very 
effectively, in the absence of parasites. 

• A "Price Volume Dealer" tries to deduce from the market 
data the current demand and supply schedules and thus 
decides whether the current price is above or below of 
the "true" value. From that data the dealer decides what 
bids and asks to quote. For this particular type of 
dealer it is also very easy to decide what volume to 
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sell, which may be a very valuable feature for a market, 
where volume is determined endogenously . 

• A "New Volume Dealer, " is a modification of the previous 
5 type. He looks on the past discounted volume only, 

without explicitly taking the price into account. If the 
New Volume Dealer observes that there were more buys in 
the past, he concludes that the price is probably above 
the true value. The simplest version of this dealer 
reacts bv increasing the ask price above the current 
market price. If the converse holds (there were more 
sells in the past) , the dealer reacts by lowering his 
bxd price. A slightly more sophisticated version of 
this dealer adjusts both bid and ask as a reaction to 
excess demand or excess supply (there is a parameter 
that allows to fine tune this dealer) . When excess de- 
mand/supply is small enough (a Volume Dealer has a 
benchmark to determine this) , the dealer concludes that 
the price is right and jumps inside the spread, 
essentially acting as a parasite. Surprisingly enough, 
even such a simple rule allows the population of the New 
Volume Dealers to track the prices really well. 

• A very important strategy for our model is the 

"Parasitic Dealer" strategy. Parasitic Dealers do not 
actively participate in the price discovery process; 
rather, they "hang back," waiting until a sufficiently 
narrow spread with sufficient volume appears on the 
board. They then "steal the deal," buying and selling 
the security inside the spread. Such strategies have 
been shown to be both effective for the dealer employing 
them, and destructive to the market's function of price 
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discovery as explained in Friedman, Daniel, and John 
Rust, eds . (1993) : The Double Auction Market: 
Institutions , Theories and Evidence, Reading, 
Massachusetts: Addison-Wesley . On at least one major 
exchange, this strategy is well understood and strongly, 
albeit unofficially, discouraged. On the AMEX, this 
process is known as "steenthing the market", named from 
the practice of shaving a sixteenth of a point off the 
existing spread to steal the deal. While it is not 
officially outlawed, it is unofficially forbidden and 
occurs very rarely in- practice. 



• A somewhat more sophisticated version of a parasitic 

dealer, called "New Parasite," also looks at the excess 
15 supply/demand and decides whether it makes sense to jump 

inside the spread or not. This version of the parasite 
performs quite well in most of the markets we have 
experimented with . 



A "Matching" dealer who contains a collection of 
observations and actions. He learns connections between 
observations and actions which are profitable. Examples, 
of observations are: "my inventory is high" , "the spread 
is large", "my bid is below the market", "my inventory 
is near zero", "trade volume is high". Examples of 
actions are more obvious: "raise my quotes" , "increase 
my bid", "narrow my spread". 

A "Classifier" dealer is similar to the Matching dealer, 
except he learns over patterns of observations. For 
instance, instead of just matching one of the 
aforementioned observations. The classifier dealer might 
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learn that "when volume is high" and "my inventory is 
near zero" and "my spread is large" it is profitable to 
"narrow my spread" . This space is much richer, clearly, 
but therefore much larger and harder to learn in. 



A "Dynamical System" dealer uses a discrete dynamical 
system to set his bid-ask spread and mean price. If one 
regards mean price, P. as a source or a sink of 
inventory, I , and the bid-ask spread, S, as a source or 
a sink of cash flow, C, then generically a system of the 
form : 



dS ( dC^ 
g 



dt 



V dt J 



20 



where f is monotonically increasing in I and g is 
increasing in C, will yield dynamics that are similar to 
those used by a real dealer. This can be seen in the P 
equation by considering the effect of changing inventory. If 
f(I 0 ) = 0, then P will remain fixed; as I increases past I 0 , 
25 will start to decrease, which will generically lead to more 
sells than buys, reducing inventory until the equilibrium 
value I 0 is again reached. Similarly, as cash flow increases, 
this implies that there are additional profit opportunities, 

(d£\ 

so the dealer should increase his spread; S\ j t j 



P 



will also 



30 

increase, causing the dealer to raise spread to capitalize on 
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the increasing market volume. Clearly, one can make this 
simple formula more sophisticated by using a predictive model 
of C or I based on the order flow and pricing observed in the 
real market . 

5 The names of these strategies are somewhat 

arbitrary and even misleading: for instance, the "New Volume" 
dealer turned out to be a very effective parasite. 

There are two senses in which evolutionary agents 
have been investigated in this model. First, in the simpler 
10 sense, running the simulation for long periods of time, while 
incorporating mutation of the parameters inside dealers' 
strategies, allows one to discover good values of those 
parameters (although the actual good values of these 
parameters may depend upon the other dealers in the 
15 simulation) . 

Second, one can use evolutionary approaches at a 
higher level to discover good strategies. Two of our dealer 
strategies are of this form: the "Matching Dealer" and the 
"Classifier Dealer." These are significantly more 
20 sophisticated types of strategies: the dealer has a whole 
repertoire of possible strategy building blocks, and by 
trying them out attempts to learn what the optimal strategies 
are. Our "Classifier Dealer" uses an approach very similar 
to classifier systems in the machine learning literature. 
25 These strategies must- be run for very long periods of time to 
allow the agents to learn the connections between observables 
(quote changes, trade volume, inventory level,...) and useful 
responses (changes to bid and offer) in the market. 

Our results indicate that the choice of strategy 
30 building blocks is of vital importance in allowing 

evolutionary learning approaches to function effectively 
(analogous results are well known in the fitness landscape 
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and optimization literature) . The literature supports our 
result that a sophisticated strategy is by no means certain 
to be more successful. For instance, the Double Auction 
Market explained in Rust, John, John Miner, and Richard 
5 Palmer (1993) : "Behavior of Trading Automata in a 

computerized Double Auction Market" found that the simplest 
parasitic strategy was the most effective. Furthermore, as 
suggested by Darley, Vincent M . , and Stuart A. Kauffman 
(1997) "Natural Rationality" in The Economy as an Evolving, 
10 Complex System II, ed.by W. Brian Arthur, David Lane and 
Steven N . Durlauf, Reading, Massachussetts : Addison- 
.Wesley (the contents of which are herein incorporated by 
reference) , there may be some bound on the sophistication of 
a strategy for it to be effective in a system of co-evolving 
15 agents. We will discuss the issue of more complex, 
sophisticated Dealer strategies later. 

Our primary results concern the influence of tick 
size and parasitic strategies on the market's effectiveness 
at price discovery. Within the current model, the underlying 
20 value of the single stock is exogenously specified, which 
allows a precise measurement of how effectively the market 
follows that underlying value, as we will discuss later. 

Our particular focus has been to investigate the 
effects of introducing parasitic dealers into the market, to 
25 determine what level of parasitism the market can support 
before its functionality becomes degraded. We then 
investigate how that degradation occurs at varying tick 
sizes . 

Our first observations are for systems consisting 
30 entirely of non-parasitic strategies . It has been discovered 
in the simulations that these systems track the underlying 
value somewhat more closely with larger tick sizes. FIG. 3 
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shows the behavior of the population of Basic Dealers for 
different tick sizes. Visually, we can see that the three 
curves plotted (best bid, best offer, and true value) are 
closer together at a larger tick size, although the 
5 difference is small (more precise quantitative measurements 
are given later) . 

The second set of observations concerns the 
introduction of parasites into the system. A number of 
different parasitic dealer strategies have been investigated, 
10 of varying sophistication. These results hold for each 

parasitic type. We observe that price discovery is largely 
unimpeded at large tick sizes, but dramatically affected for 
small tick sizes. FIG. 4 displays the results of the dealer 
mediated securities market model for a market consisting of 7 
15 Basic Dealers and one parasite. 

FIG. 5 displays the results of the dealer mediated 
securities market model for a market with a mixed population 
of dealers (we will call it a mixed market in the future) ; 
some parasitic and some not. This market contains a much 
2 0 higher proportion of parasites than the preceding ones, 
including some of the more extreme types of parasite. 

Since we are using a "true" value for the price, we 
are able to quantify the impact of different strategies and 
market structures. In particular we have examined data on 
25 the statistical distribution of (best bid - true value) 2 , 

(best offer - true value) 2 , average value of the spread and 
standard deviation of the spread over time as estimates of 
the market's ability to track price. 

The statements about price discovery and tracking 
30 the market in the preceding discussion are intuitively clear. 
A quantitative criterion, nevertheless, would be useful, 
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especially to distinguish similar cases. We will define mar- 
ket tracking as consisting of two different components: 

1) The market bid and ask prices are close to the true 
value, and 

5 2) The true value is inside the market bid and ask. 

As one can see, these two conditions do differ: it is 
possible for the market to "follow" the true price while the 
true price is outside the spread. It is also possible to 
observe that the market reacts very slowly to the true price 

10 although the true price stays inside. the bid-ask spread, 
essentially because the spread size is large. . Obviously, 
some regimes of tracking the price are more beneficial for 
Dealers and others for Investors. 

With this in mind, we quantified the market's 

15 ability to track prices by four different parameters : 

1) the square root of the average squared (bid - True 
Value) . We will denote it RMS (bid-True Value) 
where RMS means Root -Mean- Square . 

2) the square root of the average squared (ask - True 
20 Value), which we will denote as (RMS (ask- True 

Value) ) 

3) average value of the spread 

4) standard deviation of the spread. 

While a combination of these variables is not necessarily the 
25 best measure of market success (for one thing, it does not 
gives us a proportion of time when the true value was inside 
of the spread) , it gives us a very good idea of market 
performance . 

FIG. 6 shows graphs of RMS (bid - True Value) and 
3 q RMS (ask - True Value) as a function of the tick size for 

different populations of Market Makers in the dealer-mediated 
securities market model. It can be clearly seen that both 
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variables generally increase with the addition of the 
parasites to the market, and decrease with the tick size. 

FIG. 7 shows the average value of the spread as a 
function of the tick size for different combinations of the 
5 market makers in the dealer-mediated securities market model. 
There is one simulation where the spread was negative; that 
happened probably because of a very large tick size and of 
the fact that the market makers make their decisions 
simultaneously, therefore there is a possibility to 

10 "overshoot " . 

FIG. 8 shows a graph of the standard deviation of 
the spread size as a function of the tick size for different 
populations of Market Makers in the dealer-mediated 
securities market model. The results are mixed, but quite 

15 surprisingly one can see that the standard deviation 

increases somewhat with the tick size. We can also see that 
the introduction of parasites into the market can increase 
the standard deviation of the spread. 

The model is reasonably complex, and a number of 

20 other interesting phenomena have been noted (like patterns of 
wealth distribution among dealers, different market regimes). 
Our ability to measure how well the market is executing the 
process of price discovery under different circumstances can 
also help us to understand the interaction of inventory, 

25 liquidity (implied by price-volume combinations listed on the 
board) , and volatility. 

Further, we note that when investors are largely 
uninformed (i.e., their knowledge of the underlying value of 
the stock is very poor) , any collection of dealers has great 

30 difficulty in tracking the true price. Empirically, the 

reason for this is clear: when most investors are uninformed, 
no matter where the market bid and offer are located, roughly 
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equal numbers of investors will buy and sell, so very little 
can be learned by the dealers. 

One relatively obvious way in which parasitism 
affects market behavior is that, by their mere existence, 
5 parasites decrease incentives for the market participants to 
aid in the process of price discovery (Rust, John, John 
Miner, and Richard Palmer (1993) : "Behavior of Trading 
Automata in a computerized Double Auction Market", in The 
Double Auction Market) . A less obvious observation, however, 

10 is that parasites can produce useful information just by 

jumping inside the spread. The quantity and quality of this 
information does depend on the tick size. For example, if a 
parasite has undercut the market by a significant amount, the 
effects of that undercut (for instance, on market volume on 

15 the buy or sell side) tell the market participants much more 
about the liquidity, demand, underlying fundamentals, etc., 
than when the price increment is negligible. Undercutting by 
a significant amount also carries more risk and less profit. 
It is this effect, we hypothesize, which allows price 

20 discovery to continue reasonably well at a large tick size, 
even in the presence of large numbers of parasites: by 
jumping inside the spread by a large amount, a parasite 
causes a significant change in the volume of trades on that 
side of the spread, which provides useful information to the 

25 market. When tick-size is small, the marginal impact of 

jumping inside the spread on trade volume is negligible, and 
hence parasites operate without providing useful information 
to the market. A decreasing tick size will make parasitic 
behavior more profitable, since the points shaved represent 

3 q higher profit for the parasite. 

If changes in the underlying market are quite large 
(i.e., the volatility of the security is large), it may be 
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difficult for the market to track price effectively. In this 
case, it may be true that a larger tick size encourages 
better price discovery because the coarse-graining effect of 
the tick size allows for more rapid price tracking. One 
5 would anticipate that this result holds true at least for 
simple strategies, but that more sophisticated strategies 
could perhaps follow price changes equally well with a large 
or a small tick size. Our results so far support the idea 
that simple strategies can follow price slightly more 

10 effectively with a large tick size. They also suggest that, 
even if sophisticated strategies could track price well at 
all tick sizes, in the presence of parasites this tracking 
would be highly impaired at a small tick size. 

We further conjecture that a somewhat unexpected 

15 relationship between the tick size and the spread may hold if 
the ability of the system (and therefore of an individual 
dealer or investor) to adjust depends on the tick size. 
Under that condition, the dealers' profit (which depends upon 
their ability to adjust) will depend upon the tick size; 

20 since dealers will often adjust their spread to account for 
profit/risk tradeoffs/ their spread will in turn depend upon 
the tick size. 

The present invention includes an Agent -Based Model 
of a stock market with a reasonably diverse population of 

25 dealer strategies which has allowed us to investigate the 
focal questions of tick size and parasitism quite 
effectively. Our main findings at the present stage pertain 
to the effects of tick size and of parasitism on the market's 
performance and on its ability to carry on the process of 

3 0 price discovery. It has been shown in the simulations that 
1. the ability of the market to track price increases 
somewhat with tick size, and 
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2. parasitic strategies impede the process of price 
discovery, especially for a small tick size. 
While unexpected and somewhat counterintuitive, these results 
seem to be very robust in the present model . They have been 
5 repeated time and again for different combinations of 
Dealers, different levels of Investor informedness, and 
different parameters of price fluctuations. 

The market structure further comprises: a public 
limit order board must be included, together with a simulated 
10 electronic execution system; transaction costs must be 
included; and the typical statistical volume and size of 
trade must be configurable (to represent different types of 
Nasdaq stocks). This in turn will require a more complex 
range of investor and dealer strategies. First, investors of 
15 varying types should be modeled: from small consumers to 

institutional investors, and of course day- traders (who would 
probably represent the main parasitic source in the more 
realistic market) . Here an interesting addition to the 
concept of investor "informedness" would be information time 
20 lags: some investors are watching trading screens; others 

just read the Wall Street Journal. Second, dealer strategies 
account for the larger market environment (limit orders, 
greater range of trade size, etc.). Other possibilities which 
could be incorporated in the model are to allow the investors 
25 to choose between the single stock and a surrogate of the 
market index. 

There are two basic modes of operation for the 
simulation: graphical and batch. The batch mode is most 
useful for data-gathering purposes. It runs the simulation 
30 entirely without a graphical user interface, outputting to a 
file whichever data the user requests. These files are in a 
standard text format and can be read by any data/graphing 
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application for analysis and plotting. The batch version 
runs several orders of magnitude faster than the graphical 
version, and hence is ideal for the data-gathering purposes 
necessary to obtain statistically significant, robust 
5 experimental results. 

This section will first explain the control of the 
simulation via the graphical user interface, followed by the 
batch mode features. 

The graphical version of the simulation is run by 

10 double-clicking on the "BiosMarketll " icon. This will bring 
up the main simulation control window. FIG. 9 displays the 
main simulation control window of the dealer-mediated 
securities market model. 

The main window is split into two main areas: the 

15 "Quote Montage" on the left, and a recent historical view of 
the best quotes and trades as they occur on the right The 
quote montage contains each dealer's bid and each dealer's 
offer, ranked vertically from best to worst (so the largest 
bid and smallest offer are at the top of the montage) . 

20 Dealers are identified by a name (a three-letter acronym) and 
a unique color. The graphs of the dealer's cash, assets, or 
inventory (discussed later) are also color-coded for easy 
identification. The first two letters of the acronym 
represent the dealer's strategy type (for example x BD ' for 

25 Basic Dealer, X MV' for New Volume dealer, *PD' for Parasitic 
Dealer, . . . ) 

The plot on the right contains, by default, a trace 
in pink of the best offer at each time-step, a trace in yel- 
low of the best bid at each time-step, and asterisks in green 
30 of the trades that occur. Since, within each individual 

time- step, the available best bid/offer may change as trades 
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come in, not all of the trades plotted will lie directly on 
either the bid or the offer traces. 

There are five other parts to the main window: the 
current time indicator in the top left (since the simulation 
5 does not consider the opening-closing trading cycle, the 
units of time are somewhat arbitrary); the four control 
buttons in the bottom right; a display of market observables 
in the bottom left; the menus across the top of the window; 
and the ticker tape across the bottom that records recent 
10 trades. 

The four control buttons are as follows: "Quit" 
terminates the simulation; "Step" runs the simulation for a 
single time-step; "Run" runs the simulation continuously; 
"Stop" pauses the simulation (pressing either "Step" or "Run" 

15 will continue from that point) . 

The main simulation control panel contains four 
menus. First, the "File" menu currently contains a single 
item, "Exit, 11 which terminates the simulation. FIG. 10 
displays the Options menu of the dealer-mediated securities 

20 market model. The "Options" menu has four choices. FIG. 11 
displays the Set Human Dealer Prices panel and the Set Tick 
Size panel of the dealer-mediated securities market model. 
"Set Human Dealer Prices" will, if the current market 
contains a "human dealer , " open a control window that allows 

25 the user to control that dealer and thus join in the 

simulation. By clicking with the mouse, the user can adjust 
either bid, offer, or both up or down by a single tick. 
Clicking multiple times will adjust by progressively large 
amounts. To run a simulation with a human dealer, double 

30 click on the x HumanDealerMarket ' icon. Due to the memory 
consuming nature of Java interfaces, the simulations will 
tend to run faster if only one is running at a time, so we 
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recommend quitting the current simulation before starting a 
new one . 

"Set Tick Size" brings up a small window which 
allows the user to adjust the tick size of the current simu- 
5 lation on the fly. It is used by moving a slider to the 

desired position (from $1 on the left, to $ 1/100 = 1 cent on 
the right) , and then pressing the "Save" button. 

FIG. 12 displays the "Show True Value" panel and 
the "trade history" panel of the dealer-mediated securities 

10 market model. "Show True Value" toggles the display of a 
trace showing the true underlying value of the stock in the 
main display window. By default the display of the true 
value is off, so that a "human dealer" is not given 
information to which he is not entitled. Selecting "Trade 

15 History" from the menu brings up a scrolling window of recent 
trades. Each line in the window contains, in order, the 
time, investor, action that took place, volume, dealer 
involved, and price of the transaction. Each investor is 
given a three-letter lowercase name. Each dealer is in 

20 uppercase. We can see here that dealers "NVT" and "BDF" are 
dominating the market A characteristic of our markets is that 
often a small number of dealers provide most of the liquidity 
(much as in the real market) . In this case, "NVT" is in fact 
using a parasitic strategy (of the "New Volume" variety) . 

25 FIG. 13 displays the "graphs menu" of the dealer- 

mediated securities market model. In addition to the main 
trade/bid/offer display, the simulation allows the user to 
graph a number of other observables . There are six other 
built-in graphs. In addition, the simulation can be 

3 0 configured to save data to a file. This data can later be 
plotted or analyzed using any graphing/statistics package. 
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FIG. 14 displays the "cash bar graph" panel of the dealer 
mediated securities market model. The "Cash Bar Graph" 
displays a histogram of each dealer's current cash position. 
The color of each column identifies the dealer, and the 
5 height shows the current amount of cash the dealer has. 

FIG. 15 displays the "Inventory Bar Graph" panel 
of the dealer-mediated securities market' model. The 
"Inventory Bar Graph" displays a histogram of each dealer's 
current inventory position. The color of each column 

10 identifies the dealer, and the height shows the current 

volume of stock the dealer has in his possession (which could 
be negative if the dealer is short in the market) . 

FIG. 16 displays the "Supply and Demand Graph" 
panel of the dealer-mediated securities market model . The 

15 "Supply and Demand Graph" displays the supply and demand 
curves implied by the dealers' current quotes. The true 
value of the stock is also displayed. If the dealers are 
tracking the true value effectively, the graph should look 
something like the one in FIG. 16. 

2 0 FIG. 17 displays the "All Dealer Cash" graph of the 

dealer-mediated securities market model . The "All Dealer 
Cash" graph contains a color- coded trace for each dealer's 
current cash holding. It allows the user to compare the cash 
positions of the different dealers over a recent historical 

25 period. FIG. 18 displays the "Dealer Assets" panel of the 
dealer-mediated securities market model. The "All Dealer 
Assets" graph contains a color-coded trace for each dealer's 
current assets. It allows the user to compare the total 
assets of the different dealers over a recent historical 

3q period. The assets are calculated by adding the dealer's 
cash position to the value of its inventory. This value is 
calculated using the instantaneous true value of the 
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underlying stock, and hence fluctuates as that stock value 
fluctuates. The dealers themselves are of course not privy 
to this information . 

FIG. 19 shows the Dealer Inventory graph of the 
5 dealer-mediated securities market model. The "All Dealer 
Inventory" graph contains a color-coded trace for each 
dealer's current inventory position. It allows the user to 
compare the inventory of the different dealers over a recent 
historical period. In the example in FIG. 19, the " NVG " 

10 dealer went very short in the market, but then bought 

significant amounts of stock in the 84-10 0 time frame to 
balance his books , 

Thi s section describes a simple example run that 
illustrates how the market ' s price discovery process fails in 

15 the presence of parasitic dealers when the tick-size is 
small . To start the simulation, double- click on the icon 
"BiosMarketExample . " After a brief delay, the simulation 
interface appears. Before running the simulation, to see the 
process of price discovery in action, it is best to display 

20 tne underlying value in the bid/of fer/ trade graph to the 
right of the main control panel. To do this, go to the 
"Options" menu and select "Show true value." To see the 
results most dramatically, it is best to change the tick-size 
first to a very large value, and then to a very small value. 

25 So, on the "Options" menu, select "Set Tick Size" and click 
on the blue slider to move it to the extreme left, so that 
the text entry reads "1." Now click on "Save." The tick-size 
in the market will now be set to $1; i.e., a very large 
value . 

3 Q Now click on the "Run" button in the main control 

window, and watch the simulation unfold. This simulation 
contains a single parasitic dealer, of the "New Volume" 
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variety . FIG . 2 0 displays a panel showing how the dealers 
track the true value reasonably well in the dealer -mediated 
securities market model. 

Once you are content that price discovery is 
5 functioning well, set the tick-size to a very small value. 
There is no need to stop the simulation to do this; simply go 
back to the "Set Tick Size" window and move the slider to its 
rightmost end so that the text entry reads "100" (or 
something near 100). Again, you need to click on the "Save" 

10 button for the change to take effect. FIG. 21 displays a 

panel showing the simulation results for a smaller tick size 
in the dealer-mediated securities market model. Price 
discovery is not as effective . With larger numbers of 
parasites, or a single more extreme parasite, the tracking of 

15 true value fails almost completely. 

Almost all of the data we examined to draw our 
conclusions were generated by running the simulation in batch 
mode . This has several advantages : first and foremost , the 
simulation runs several orders of magnitude faster; it is 

20 also possible to pre-specify a total number of time-steps for 
which to run it (anywhere from 5000 to 500, 000 were used in 
practice) , and it is easier to run several simulations in 
parallel . 

It is possible, to run the simulation in these two 
25 different ways because the core of the simulation is highly 
modular. It is necessary simply to connect that core with 
either the GUI or the batch code to run the simulation. To 
run a simulation in batch mode, one simply types (from the 
command line) : 

30 

java BatchMarket [prefix] [duration] 
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where " [prefix] " is a character prefix the simulation uses to 
identify input and output files for this run, and 
" [duration] " is the number of time-steps for which to run the 
simulation (if omitted, this defaults to 10,000 time-steps). 
5 For instance, to run a simulation without the GUI for 5000 
time -steps , using the same input files as the example run 
described earlier, one types: 

java BatchMarket EG 5000 

10 

The files that describe the initial conditions (dealer mix, 
etc.) for the example run are called " EGdealerMix . dat , " 
"EGoutputParms.dat, " and "EGplayerParins.dat." 

The user can configure the initial conditions for 
15 the simulation by editing those files. The " output Parms " 
file contains "name value" pairs in the following format: 

outputDealerCashData true 

output lnve s tor Invent oryDat a false 

20 

The possible parameters to configure are as follows: 

• base_f ilename: the suffix to use for all data files 
created by this simulation run 

25 • outputData: true/ false flag that globally turns data 

output on or off. If off, it overrides all other data 
output flags . 

• OutputlnvestorCashData: if true, investor cash holdings 
are output to file at each time- step 

30 • OutputDealerCashData: if true, dealer cash holdings are 
output to file at each time- step 
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• OutputlnvestorlnventoryData : if true , investor inventory 
holdings are output to file at each time - step 

• OutputDealerlnventoryData : if true , dealer inventory 
holdings are output to file at each time-step 

5 • OutputValueDebugData : if true, best bid, best offer, 

and true value are output whenever a trade takes place. 

• OutputObjectFilesRate: rate at which to save (using 
Java's serialization capabilities) all Dealer objects to 
file for later examination 

10 • OutputAssets: if true, dealer assets are output to file 
at each time-step 

• OutputBids: if true, each dealer's bid is output to file 
at each time- step 

• OutputOf f ers : if true , each dealer ' s offer is output to 
15 file at each time-step 

• OutputBankruptcy : if true, data on which dealers and 
Investors go bankrupt is output to file 

• Unif ormRandomWalk: if true, using a uniform random walk 
with drift, otherwise a log-normal random walk 

20 • FoolishnessMax: the standard deviation of the investors' 
knowledge-error in the true value 

• DisallowArbitrage : if true, then disallow dealer quotes 
that would allow obvious arbitrage opportunities 

• RandomWalkSeed: can be used to fix the seed for the 
2 5 random number generator used by the random walk 

• TickSize: sets the initial tick size. 



The "dealerMix" file contains the list of dealers which take 
part in the simulation. Usually we fix this at eight 
3 0 dealers. The file simply contains a list of the names of the 
Java classes which implement the different dealer types. For 
example " EGdealerMix . dat " contains the following list : 
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BasicDealer 
BasicDealer 
BasicDealer 
BasicDealer 
5 BasicDealer 
BasicDealer 
NewVolumeDealer 
NewVolumeDealer 

10 The third file, "playerParms " contains default initial 

parameters (including the inventory and cash holdings) for 
the different dealer types. It contains a number of lines, 
the first item on the line being the name of the dealer type, 
followed by the number of arguments required by the object's 

15 constructor, followed by the arguments themselves. For. 
example, the line for the basic Parasitic Dealer is: 

ParasiteDealer 5 1000. Of 100 3 2. Of 3 

20 So the dealer is initialized with $1000 in cash, 100 units of 
stock, a responsiveness of 3 (the scheduler will allow the 
dealer to act every three time steps on average) , a risk 
aversion of 2.0 and will accept a maximum of 3 trades inside 
the spread before moving his quotes outside. 

25 This aspect of the present invention has been 

developed using an object-oriented simulation of a single- 
security market with the Java programming language. Three 
principal classes make up the simulation: the Market class, 
the Dealer class, and the Investor class. In addition, there 

30 is a support class, the GuiMarket class, which provides a 
user interface and ties the other objects together in an 
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actual simulation . Individual instances of these classes are 
obj ects . 

• Market: The Market class provides a venue where the 

5 Dealer and Investor agents can interact. The Market 

provides the following functionality: 

• A time-step method, which runs the simulation 
forward one step in time . 

• Methods for matching Investors who wish to buy or 
10 sell the security with the Dealers who are trading 

the security at the Investors' desired price. 

• A Quote Montage object, which displays the current 
bid and ask prices and volumes of all Dealers in 
the market. It also provides a method for 

15 displaying the current best bid and ask and their 

volumes . 

• A Scheduler object. The scheduler determines the 
order and frequency of trades made by each Dealer 
and Investor. It allows each player waiting on a 

2 0 queue to trade, and then puts them back onto the 

queue to wait for their next opportunity. 

• Dealers and Investors: These are both subclasses of a 
Player superclass. The Player class provides interfaces 

25 for keeping track of inventory and capital and executing 

trades. Dealers have the additional ability to interact 
with the Market's Quote Montage object, while Investors 
are able to query the Market for the current "true" 
value of the security. Both Dealers and Investors have 

3Q strategy functions that determine their trading desires. 

Investors decide whether to buy or sell the security, 
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while Dealers decide what bids and asks to post on the 
Market's Quote Montage. 

• GuiMarket : The GuiMarket provides the user interface, 
5 and the GuiMarket object is the entry point into the 

simulation when run in graphical mode. The main 
functions of the GuiMarket object are: 

• Instantiation of the Market object with desired 
initial conditions . 

iO • Display of the Market's Quote Montage object as 

trades occur. This panel displays pertinent in- 
formation about the time and spread. It also 
displays all the dealers, the volume, and the bid 
price on the left side , and the dealers , volume, 

15 and price offered on the right side. 

• Display of Buy/Sell graph as trades occur. The 
volume of stock for current bid prices is graphed 
against the volume of stock at the current offer 
price. 

20 • Buttons for running/stopping/stepping/quitting the 

simulation. 

• Instantiation of a Ticker object, which currently 
displays to the scrolling ticker tape at the bottom 
of the main control window. The ticker records 

25 market activity, i.e., time, Investor name, buy or 

sell, volume, dealer name, and price. 

• Initialization of a Config object 

• Config : this obj ect provides the following 
30 functionality: 

• Market parameter set-up, which includes 
initialization of the number and types of dealers 
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(BD - basic dealer, PD - parasitic dealer, etc) and 
number and types of investors (Fl - foolish 
investor, etc.); setting the number of time 
iterations for a simulation, setting initial cash 
5 and inventory for investors and dealers and their 

response time statistics; setting the "true" price 
of the stock. These parameters are accessed 
through the three parameter files discussed above. 

10 • BatchMarket : This provides the non-graphical entry point 
to the simulation. It initializes a Config object which 
reads in the input files (simulation run parameters) , 
sets the market running for the specified number of 
steps, and then exits. The input files specify which 

is data is to be saved during the run. 

When the market is run (either via M java GuiMarket" or "java 
BatchMarket") another command line parameter is required. 
This parameter, as described previously, is the prefix used 
20 for the three input files. So, to create a new simulation, 
one simply edits those files, and then executes the command 
"java GuiMarket <prefix>" or x java BatchMarket <prefix> 
<duration>" . 

FIG. 22 displays a panel of results obtained by 
25 executing the dealer mediated securities market model. 

Predicting Market Volatility 

The space of economic models 10 5 further include 
models which are particularly adapted to predict volatility. 
3 q It is clear that current efforts to predict market volatility 
based on the Black- Scholes formula for option pricing do not 
work particularly well. 
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As is well-known the implied volatility obtained by 
running the formula backwards, using stock and option prices, 
varies nonlinearly with the difference between the current 
price of the underlying and its exercise price, instead of 
5 being independent, so there is something wrong with the 
structure of the Black-Scholes model. 

It is also well-known (and possibly related) that 
parameterising market returns by the lognormal density does 
not work very well for short term data, i.e. the density 

10 p(s , /s) where s is the asset value now and s 1 its 

value in the future, is not lognormally distributed if the 
future is not more than about two weeks, but follows a 
Pareto-Levy distribution. 

It follows that better formulas are need to predict 

15 volatility both within the framework of implied volatility, 
by developing better models of the variation of stock returns 
than the simple Markov process which underlies the 
Black-Scholes equation, but also by developing better ways to 
analyze the data provided by market time series . 

2 o We propose to start the process of finding better 

market volatility predictors by constructing a nonlinear 
state space model (NLSSM) of the market. Such a model 
contains hidden state variables and parameters. To 
estimate these variables and parameters we propose to use a 

25 version of the expectation-maximization or EM-algorithm. 

Such an algorithm has two steps (a) an expectation or E-step 
in which hidden variables are estimated from 

given parameters, and (b) a maximization or M-step in which 
parameter estimates are updated given hidden variables. The 
3Q E-step can be carried out using a Kalman filter, which can be 
viewed as a generalization of the backpropagation algorithm 
of neural networks . 
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It is clear by now that in applying such algorithms 
to markets, data should be as low dimensional as possible. 
There are now a number of differing techniques for dimension 
reduction which can be used. There is a class of 
5 unsupervised learning algorithms (as distinct from 

backpropagation which is a supervised learning algorithm) 
embodied in neural networks available for this purpose, e.g. 
the principal components analysis algorithm (PCA) and various 
variants (localised PCA, nonlinear PCA) , and a 

10 close relative, the independent components analysis algorithm 
(ICA) . These algorithms can be used to find optimal 
representations of the data before setting up the NLSSM. 

Several other issues associated with this approach 
suggest themselves . 

15 (a) What is the relation of the hidden variables in an NLSSM 
to those found in causal or agent based models of the market? 
Evidently better market models can give rise to better state 
space descriptors. Such descriptors could include cooperative 
interactions between traders, such as herding. 

2 0 (b) There are also many new developments associated with 

parametric approaches to market modeling that use Pareto-Levy 
distributions etc, or nonlinear Langevin models, e.g. 
Bouchoud et.al. and Farmer. These can also 
give rise to better state space representations. 

25 (c) Finally it is also possible that ways can be found to 
optimize NLSSMs using search algorithms, e.g. genetic 
algorithms (Gas) to find better hidden state variables and 
their combinations. This will require a measure of 
improvement in the ability of such models to predict market 

30 data or properties such as volatility. 

Optimization 
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Once the causal models have been evolved and 
adapted to fit real market data, they are used to generate 
data more widely for a multitude of situations, including 
what -if scenarios- In particular, they can be used: 
5 1. for predicting volatility; 

2. for evaluating and predicting risk for any desired 
portfolio, allowing ranking of different 
portfolios ; 

3. for optimization of portfolios; 

10 4. for pricing of options and other derivatives; 

5 . for choice of optimal hedging strategies . 



Because the models produce non-Gaussian 
distributions of return, a single simple measure of risk such 

15 as volatility is not adequate to characterize risk, rank 
portfolios, optimize portfolios, price options, and hedge 
risks. Practical risk measures that we use include the 
probability of loss greater than a specified threshold, the 
expected worst drawdown in a given period, the expected time 

20 between significant drawdowns, and the value at risk (VaR) at 
a given confidence level. We extract these quantities from 
our models using combination of analytic and statistical 
means . 



Minimizing Values at Risk 

The present invention includes techniques to 
minimize the value at risk of a portfolio. Value at risk is 
a single, summary, statistical measure of possible portfolio 
30 losses. Specifically, value at risk is a measure of losses 
due to "normal" market movements. Losses greater than the 
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value at risk are suffered only with a specified small 
probability . 

Using a probability of x percent and a holding 
period of t days, a portfolio's value at risk is the loss 
5 that us expected to be exceeded with a probability of only x 
percent during the next t-day holding period. 

The technique to minimize the value at risk uses 
historical simulation. Historical simulation .is a simple, 
atheoretical approach that requires relatively few 
10 assumptions about the statistical distributions of the 

underlying market factors. In essence/ the approach involves 
using hxstorical changes in market rates and prices to 
construct a distribution of potential future portfolio 
profits and losses, and then reading off the value at risk as 
15 the loss that is exceeded only x percent of the time. 

The distribution of profits and losses is 
constructed by taking the current portfolio, and subjecting 
it to the actual changes in the market factors experienced 
during each of the last N periods. That is, N sets of 
20 hypothetical market factors are constructed using their 

current values and the changes experienced during the last N 
periods. Using these hypothetical values of market factors, N 
hypothetical mark-to-market portfolio values are computed. 
From this, it is possible to compute N hypothetical mark- to- 
25 market profits and losses on the portfolio. 

The following discussion describes the technique 
for isolating low value at risk portfolios. Let us consider 
a single instrument portfolio, in this case stocks traded on 
the New York Stock Exchange and Nasdaq . market s . For this 
3 Q instrument, there exists tremendous amounts of data. If we 
assume a one day time horizon ( t= 1), then the data we are 
interested in are the daily closing prices of every publicly 
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traded stock on the two markets. Such data exists for 
thousands of stocks for tens of thousands of days. From 
these data, it is possible to construct an m x n matrix 
(where mis the number of stocks, and n is the number of days) 
5 of prices . 

Let us assume that within this collection of 
stocks, there are, pairs, triplets, quadruplets, etc., of 
stocks whose values at risk are lower as a group than any of 
the stocks individually. This occurs because sets of stocks 

0 whose price changes are anti-correlated will have low values 
at risk. When the price of one stock goes down, the price of 
the other tends to go up. The chance that both stocks go down 
together is lower than the chance that two stocks chosen at 
r~random would go down together because the stocks are anti- 

5 correlated. This reduces value at risk. 

The optimal portfolio would group ant i -correlated 
stocks in the optimal proportions to minimize value at risk. 
Because there are so many stocks, however,, the space of all 
possible portfolios is too large to search exhaustively. 

^ Gene t:ic algorithms are well suited to finding good solutions 
to just this type of problem in reasonable amounts of time. 
The algorithm works as- follows: 

Step 1 : 

Start with m portfolios. Each portfolio can be represented as 
5 a vector of length m. Each bit (m A ) in the vector is either a 
1 or a 0 signifying that the i th stock is either included or 
excluded from the portfolio. This can later be extended to 
letting each bit specify the number of shares held rather 
than simply inclusion or exclusion. To each portfolio, assign 
0 a random number of stocks to hold such that every possible 
portfolio size is covered (at least one port folio excludes 
all but one stock, at least one portfolio excludes all but 
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two stocks, ... at least one portfolio includes all the 
stocks) . Once the number of stocks to hold has been assigned, 
let each portfolio randomly pick stocks until it has reached 
its quota. 
5 Step 2 : 

Go back in time n/2 days (halfway through the data) . For each 
of the m portfolios, compute the value at risk for the nl2 
days that precede the halfway point. 
Step 3 : 

10 Randomly pair portfolios. For each pair of portfolios, let 
the portfolio with the higher value at risk copy half of the 
bits of the lower value at risk portfolio (i.e. randomly 
select half of the bits in the more successful portfolio. If 
a bit is different, the less successful portfolio changes its 
5 bit to match the more successful portfolio) . The portfolio 
with the lower value at risk remains unchanged. 
Step 4 : 

Repeat steps 2 and 3 until some threshold for value at risk 
is achieved . 

0 

In this way, clusters of ant i -correlated stocks will tend to 
spread through the population of portfolios. The hope is that 
this method will ultimately select for most or all of the 
good clusters. Notice that this method may also alight upon 
5 the optimal number of stocks to hold in a portfolio. For 

example, if the minimum value at risk portfolio contains only 
three stocks, three-stock portfolios will tend to propagate 
through the population. 

t0 Additional Techniques for the Analysis of Risk 

The present invention includes additional 
techniques for the analysis of risk. The general 
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understanding of portfolio risk requires an understanding of 
three contributing problems. The current understanding of 
these three problems is insufficient to accommodate the 
challenges posed by modem portfolios. The first problem is 
5 volatility. It has long been known that Gaussian 

approximations to volatility do not correctly describe the 
behavior of markets, and that price fluctuations show long 
tails. This means that large deviations are much more likely 
than conventional theory suggests. The second issue is that 

10 of interdependence. In many areas of interest, elements of a 
portfolio do not move independently of each other, but rather 
influence each other in ways both subtle and complex. Current 
methods only uncover a minimal rendering of this complex 
structure. The third issue is that of time dependence. Many 

15 portfolios contain elements that do not mature on the same 
time scale, but are nonetheless dependent. Again, 
conventional portfolio analysis and optimization techniques 
do not address the subtleties of interacting time scales. 

It was originally pointed out by Mandelbrot in 1961 

20 that price fluctuations in speculative markets follow a 
particular class of fat-tailed (relative to Gaussian) 
distribution known as a Levy distributions. Levy 
distributions can be created by correlated random walks, and 
have been studied in statistical physics under the rubric of 

25 anomalous diffusion. In the last two years, these 

distributions have been revisited in the area of option 
pricing, but not yet in the area of portfolio management. 
Since Levy price fluctuations do not add together in the same 
way as Gaussian fluctuations, results for expected portfolio 

3 Q risk and associated volatility will in general be different. 
The present invention incorporates an underlying Levy model 
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of volatility for portfolio risk analysis which accounts for 
the different additivity characteristic. 

In the construction of portfolios, it is 
experimentally known that the prices of certain stocks are 
5 correlated, and this correlation is typically measured using 
a covariance matrix. The covariance matrix has two implicit 
assumptions which we believe are wrong: Fluctuations in 
prices are Gaussian (see above) and correlations between 
stocks are describable with pair-wise interactions. 
10 The present invention modifies this approach in two ways: 
1. The covariance matrix requires a large amount of data 
for accurate results. Extending the covariance matrix 
method to higher order interactions (three- or four- 
point interactions) requires an exponentially increasing 
15 amount of data. We separate the analysis of 

interdependence into effect and magnitude. The effect 
aspect is obtained by encoding price fluctuations as 
( + ,-) instead of numerical values. Now we have prices 
encoded as binary strings instead of numerical 
sequences. Since the fundamental activities of a market 
are buying and selling, and their attendant effects are 
the raising and lowering of prices, we believe that the 
binary encoding is a more fundamental signature of 
market microstructure than the actual prices. The 
magnitude of gains and losses are obtained by the 
statistics of large numbers of players making "atomic 11 
buy and sell decisions . 

Once we have encoded the market dynamics of 
individual instruments as bit strings, we have 
essentially a telegraph representation of information. 
This is amenable to the tools of Information theory, a 
field developed to analyze the dynamics of information 

- 65 - 



SUBSTITUTE SHEET (RULE 26) 



WO 00/38094 



PCT7US99/30640 



transfer in telephony. Information theory allows us to 
measure correlations at arbitrary levels of 
interconnectedness , and although the data requirements 
scale exponentially as interconnectedness increases, the 
5 constant in front of the exponent is much smaller than 

the covariance case because of the binary nature of the 
data. Interconnectedness is measured by a quantity 
called mutual information, and the assumptions 
associated with it are less stringent than the 
10 assumptions required to measure covariance, and in 

particular are not dependent on the assumption of a 
normal distribution . 

2 . The present invention uses the measure of mutual 

15 information to construct a phylogeny of interdependence, 

using the technique of minimal spanning trees coupled 
with higher order information correlations to remove 
degeneracies, (multiple solutions satisfying the same 
constraints) Once we have constructed a phylogeny (a 

2 0 tinkertoy like structure showing connections of 

influence) , we can use time- structured data to obtain 
directions of influence. This directed map allows us to 
model the propagation of financial disturbance through a 
web of connections. This is an important tool for 

25 constructing a portfolio of minimum risk, because it 

decomposes portfolio risk into an ensemble of 
interconnected contributing factors which can then be 
optimized to obtain the desired results. 
Note: The connections discussed above can be 

3q endogenously or exogenously generated. If the portfolio 

in question consists of internal assets, (R&D portfolio, 
for instance) there is some control as to the degree and 
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nature of the interconnection. Hence the optimization 
procedure is somewhat different , as the interconnections 
are no longer viewed as given but are now variables over 
which we have some control . 
5 In many areas of interest, different time scales are 

an important consideration . Examples include bonds with 
different dates of maturity, and internal portfolios with 
different payback profiles. Optimizing over such portfolios 
requires understanding the spectrum of possible paths that a 
10 portfolio can take over time, and again interdependencies and 
large fluctuations make standard approximations of Gaussian 
uncertainties inaccurate . The present invention uses 
techniques borrowed from non-standard statistics (large 
deviation theory, sampling theory) and quantum field theory 
15 (path integrals) to generate a forward curve of the behavior 
of a complex portfolio over time. Not only is the end result 
important, but the shape of the curve over time is important, 
as there are many quantities of interest which are dependent 
on the local shape of the curve. 

20 

Evaluating and Minimizing Risk 

The present invention includes additional 
techniques for portfolio optimization using sampling and 
selection to evaluate and minimize risk for a portfolio of 

25 assets with uncertain returns. Consider the general setting 
in which a holder owns portfolio of assets. The assets may 
be financial instruments (such as stocks, bonds, options, or 
other derivatives) or investments in research and development 
with unknown future payoffs (e.g. the drug leads pursued by a 

3 0 pharmaceutical company) . In this setting, where the future 
rewards are uncertain, there are two important concerns of 
the holder of the portfolio. Firstly, it is important to 
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quantify the risk (the amount of money, that could be lost) 
over some time horizon. Secondly, the holder wishes to 
structure the portfolio so as to minimize the risk. In this 
document we will focus on answering these questions for 
5 portfolios of financial instruments but the ideas are more 
generally applicable . 



asset in the portfolio. If there are N assets in the 
portfolio let xftj be the N- vector representing the values at 
10 time of all components of the entire portfolio. The value of 
the entire port folio to the holder is specified as some 
function f (x) of the values of the assets. Typically, this 
function might lie a linear combination of the asset volumes 

N 

times the asset prices, i.e.f (x) =52 V '± K ± - Furthermore let 
15 i=1 

P(x',fc'|x, t) represent the probability that the asset prices 
are x* at time t' > t given that the asset prices were x at 
time t. If t indicates the present time and x represents the 
present value of the assets then the expected value of the 
20 portfolio at some time t' in the future is 



This value indicates the expected worth of the portfolio but 
does nothing to tell us what the risk is, i.e. what we might 
conceivably lose. To get at this quantity we realize that 
from P(x' , t'|x, t) we can also determine the probability 
P(v|t)that the value at. rime t is v: 



Let x ± (t) represent the value at time t of the ith 



V{ t'|x, t) 




30 



P(v\t) 




- 68 - 



SUBSTITUTE SHEET (RULE 26) 



WO 00/38094 



PCT/US99/30640 



This probability is the fundamental quantity which allows us 
to assess risk since it gives the probabilities for all 
potential outcomes. Thus for example we can say things like 
"with 95% confidence the most money that will be lost, is 
v* . n In this case v* is determined from the requirement that 
only 5% of the time will more money be lost, i.e. 

| V/ dv P(v| t) - 0.05. 

Other measures of risk are similarly based on P ( v\ t ) . 

The risk will depend sensitively on the precise 
form of P (x', t'| x, t ) . To see this, imagine that a pair of 
assets i and j are anti -correlated with each other (i.e. when 
the price x ± increases the price usually decreases) . If we 
invest equally in both stocks then the risk will be small 
since if the value of one stock goes down the other 
compensates by going up. On the other hand if the price 
movements of assets are strongly correlated then risks are 
amplified. To evaluate and manage risk it then becomes 
0 paramount to identify set of assets that are correlated/anti - 
correlated with each other. This observation forms the basis 
of traditional value at risk analyses (VaR) .in which the risk 
is assessed in terms of the covariance matrix in asset 
prices. The covariance matrix includes all the possible 
5 pairwise correlations between assets. 

While traditional VaR captures pairwise variations 
in asset prices it completely ignores higher order 
relationships between variables, e.g. when assets i and j go 
up asset k goes down. Moreover the Gaussian assumption 
0 inherent in VaR is known to be false. What is needed is a 
more general approach. The present invention describes new 
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risk management techniques which move beyond pairwise VaR. 
We propose two methods by which this can be accomplished. 

The core of the new ideas is to recognize that 
information about higher order relationships can be uncovered 
5 by looking at the VaR of subsets of assets from the 
portfolio. 

Imagine for a moment that a specific set of assets 
covaries with each other in some predictable way. Knowledge 
of this co-variation can be used to devise a risk averse 

10 combination of these particular stocks. Since the variation 
involves all four stocks it can never be determined by only 
looking at pairs' of assets. How might important, clusters of 
stocks like these be discovered? 

The first point to note is that the historical 

15 record of asset prices and portfolio values provides a 

training set from which we can discover these clusters of 
stocks. The historical record provides a data set which 
includes the true VaR because the future value of the 
portfolio is known from the historical data. Let v represent 

20 the true VaR for a particular portfolio x at a point T into 
the future. From the historical record we can form the data 
set D= fx., v.} and thus estimate the VaR for the assets in the' 
chosen portfolio, i. e. P(v|x) . If we assume that the 
stochastic process that generated D is stationary then the 

25 same relationship discovered in D will also hold in the 

future. Once the mapping from a set of stocks to a VaR has 
been determined we can search over subsets of stocks to find . 
a combination that gives particularly low VaR. 

We begin by making the simple assumption that 

30 P(v|x)= 5 (v - /i(x)), i.e., it is characterized entirely by 
its mean value /z (x) . This mean value will differ for 
different subsets of stocks. Without much additional 
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complication we could also include the variance around this 
mean and assume that fluctuations around the mean are 
Gaussian: P (v|x) = N {fi (x) , o 2 (x) ) . Obviously from the 
data D we could try to infer much more complicated 
5 relationships but for the moment we stick with simple case. 

Given that we can determine the true average VaR 
for any set of stocks we can identify those assets within a 
portfolio of N assets that form good combinations. 
Computationally the following scheme can be used to identify 

10 good subsets of assets. We assume, that the optimal subset of 
assets is of size n « AT. Starting from the original 
portfolio randomly form portfolios of half the size by 
sampling stocks (without replacement) from the entire 
portfolio. The probability that any one of these randomly 

15 generated profiles contains all n assets is approximately 
l/2 n . Thus as we examine significantly more random 
portfolios than this it is likely we will obtain at least one 
subset containing all n assets. For each of the randomly 
generated portfolios of N/2 assets we determine its VaR by 

20 calculating it from D and keep those portfolios with high 
VaR. In this way we filter out only the most promising 
portfolios, i.e. those that contain the subset that we are 
seeking. This process can then be iterated further. From 
these remaining portfolios of size N/2 we randomly generate 

25 make portfolios of half the size (no N/4) . Assuming that, at 
least one of the size N/2 portfolios contained the desired 
cluster the probability that one of the size N/4 portfolios 
contains the full subset, is again l/2 N , We can keep 
iterating this process of generating and filtering portfolios 

3q and each time come closer to good subsets of stocks. 

After m iterations of this procedure the portfolio 
size is N/2 m . Let inbe the largest value of m such that N/2 m 
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> n (i.e., the largest portfolio that contain all n assets) 
and let m = m +1. If our hypothesis is correct then we 
should see an abrupt increase in the VaR from m to m since we 
can no longer form a risk averse combination of all the n 
5 assets at m. This fact alerts us to the important fact that 
n must lie between N/2~and N/2 m . At this point we can then 
sample from the portfolio of size N/2 m to form new portfolios 
of size (N/2 B + N/2 m )/2. The extremal VaR values of these 
new portfolios will either be comparably to the N/2 m ±n which 
10 case [N/ 2* + N/2~) /2< n < N/2 m or comparable to N/2- in which 
case (N/2 B < n < (N/2 m + N/2~) /2 . Iterating this procedure we 
can determine the optimal subset size n. Knowing the optimal 
n we can look at many different subsets of this size to 
eventually pick out the precise combination of the n assets. 
15 FIG. 23 provides a flow diagram of a method for 

evaluating and minimizing risk. In step 23 05, the method 
selects N assets and initializes the portfolio size. In step 
2310, the method generates portfolios of the specified size 
from the selected subsets. On the first iteration of the 
20 loop, the selected subset will be the set of N assets. In 
step 2315, the method computes the value at risk for each 
portfolio generated in step 2310. In step 2320, the method 
selects the optimal portfolios (subsets) from those generated 
by step 2310. 

25 in step 2325, the method determine whether any 

portfolio contains the desired properties. If at least one 
portfolio contains the desired properties, control proceeds 
to step 2330. In step 233 0, the method decreases the subset 
size before generating new portfolios in step 2310 in the 

30 next iteration of the loop. If the method determines in step 
2325 that no portfolio contains the desired properties, 
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control proceeds to step 2335. In step 2335, the method 
iteratively searches for an optimal subset size. 

There are many variations to this basic method that 
might improve its efficiency. It may not be optimal to halve 
5 the portfolio size at each step since we might want to have a 
higher probability of retaining the subset intact. The best 
number of random portfolios to generate and test can also be 
fine-tuned to make the search more efficient. Simple 
analytical model can be built to optimize these algorithm 

10 parameters. 

Another important extension is to realize that the 
sample and select method outlined above can be used to 
determine subsets with any desired properties. Nothing in 
the above procedure restricts its use to minimize VaR. The 

15 method can be used to extremize any figure of merit. Along 
these lines what happens in practice is that there are 
actually more than one objective in designing a portfolio: 
we want to minimize risk but also maximize a profit. Is 
there a way of modifying the above procedure to account for 

20 more than a single objective? At the very least we must be 
able to balance risk/reward the way it is conventionally 
done . 

There is however a very simple extension of the 
above method to handle the case of multiple objectives. Sub- 

25 sampled portfolios are generated as usually but the selection 
criteria amongst portfolios is modified. Now Instead of 
picking sub- sampled portfolios which have the best VaRs we 
measure a number of objectives for each of the particular 
sub-sampled portfolio and keep those sub-sampled portfolios 

30 which Pareto dominate all other portfolios (generated at the 
present iteration or all previously generated portfolios) . 
With this simple change in the selection criteria everything 
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else proceeds the same way. At the end of the process we 
obtain a subset which is hopefully good on all objectives. 

There are also a host of issues involving sampling 
to form the new portfolios. In the above I have assumed that 
5 the sampling is uniform. Of course this need not be the 

case. The distribution can be altered over time (much as in 
boosting) based on the VaR (or distributions over VaR) over 
previous samples. This also ties into pooling in that the 
samples may be constructed in a much more deterministic 
10 fashion. 

There are also a number of possible extensions 
based upon the fact that the mapping from the asset subset 
can be determined. At this point it remains unclear how to 
exploit this fact. The basic procedure might also have to be 
15 modified slightly to fit the context. 

One of the main difficulties of the procedure may 
occur in the initial phases of the algorithm. Initially we 
must filter out subsets that contain the desired small 
subset. There arises the signal/noise ratio issue on whether 
20 the small optimal subset can be detected in the larger 

portfolio. An estimate of this would be useful to know how 
many subset samples to generate. This has been tested and 
the results will be discussed below. 

In order to test the algorithm, a random number 
25 generator was used to create a "toy" set of stock data in 
which all of the clusters were known a priori. 

FIG. 24 provides the results of the method for 
evaluating and minimizing risk from executing on 500,000 
random draws from this "toy" world. The world consists of 32 
3 q stocks in which an investor may invest. The solid square 
plot is a histogram of VaRs for portfolios which include 
exactly 16 stocks. The mean VaR for these portfolios is 
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-0.96 (where negative means the investor loses money). The 
smaller histogram of hollow circles shows the VaRs only for 
those portfolios which, through random sampling, failed to 
include any of the good clusters. The mean for the no- 
stocks- in-clusters (SIC = 0) portfolios is -1.08, clearly 
worse than for the population as a whole. This is exactly 
the signal we need to separate the good portfolios from the 
SiC=0 portfolios. 



are selected and 8 -stock children are made from them. All 
16-stock portfolios which lie outside the 99% confidence 
interval of the SIC=0 portfolio distribution are selected. 
The number of children to draw from each is determined by the 
histograms shown in Figure 25. 



First, the total number of portfolios that contain zero 
stocks in clusters has increased dramatically. This is 
because as you remove stocks, your chance of breaking up 
clusters increases exponentially. Second, the entire 
distribution has shifted to the left, meaning that the VaRs 
have gotten worse (mean -1.66). This is because their is 
value in diversity and diversity must decrease as portfolio 
size goes down. Third, the distribution of SIC=0 portfolios 
is still worse (mean - -1.72) than the distribution of all 
portfolios . 



the number of child portfolios necessary to ensure that some 
children still contain intact clusters. The number of 
children containing clusters is negatively binomially 



At the next iteration, the best 16-stock portfolios 



There are three features worth noticing in FIG. 25. 



This third feature allows the modeler to determine 




- 75 - 



SUBSTITUTE SHEET (RULE 26) 



WO 00/38094 



PCT/US99/30640 



where n is the number of children, r is the number of 
children that contain clusters, and p is the probability of 
randomly selecting a portfolio which is to the right of the 
99% confidence interval of the SIC=0 histogram. 
5 Using these facts, a preliminary model has been 

calibrated and successfully finds clusters on more than 95% 
of the runs . 

The present invention also includes a method for 
portfolio modification. There are other methods to try to 

10 identify beneficial changes to a portfolio. Traditional VaR 
theory measures the effects of modifying (i.e. increasing or 
decreasing the holding) a position in one of the assets. As 
we have seen, if higher order combinations of assets are 
important then the effects of a single asset might be minor 

15 There is an important practical reason why traditional VaR 
focuses on the changes of only a single asset. If the 
portfolio is of size N and we consider changes involving m 
assets then on the order of N m stocks must be examines. 
Consequently, for practical reasons attention is usually 

20 restricted to m=l or single asset changes. 

In contrast, the present method determines the 
optimal number of assets to change while searching for an 
optimal portfolio. FIG. 26 displays a flowchart illustrating 
the method for portfolio modification. In step 2610, the 

25 method infers a fitness landscape from historical data. The 
landscape can be inferred from historical data using 
techniques described in the co-pending application titled, 
"An Adaptive and Reliable System and Method for Operations 
Management", U.S. Application No. 09/345,441, filed July 1, 
30 1999, International Application No. PCT/US99/15096 , filed 
July 2, 1999, the contents of which are herein incorporated 
by reference. 
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In step 2620, the method determines the optimal 
searching distance as described in co-pending international 
application No. PCT/US99/ 19916 , titled, "A Method for Optimal 
Search on a Technology Landscape", filed August 31, 1999, the 
5 contents of which are herein incorporated by reference. Once 
m is known, the method searches for optimal portfolios at the 
optimal searching distance in step 2630. For example, the 
method shown in FIG. 2 3 could be used to find an optimal 
portfolio after the method of FIG. 26 determines the optimal 

10 searching distance. 

In this case, the landscape is defined over 
portfolios. The fitness is the VaR measurement of the 
portfolio and two portfolios are neighbors if they differ in 
the holding of a single asset. So for example a portfolio of 

15 five assets might be represented as [2, 1, 1, 3, 1] meaning 
that the portfolio consists of two units of asset one, one 
unit of assets two , thee and five, and three units of asset 
four. 

20 Portfolio Management - Optimization of Risk-Revenue 

The present invention further includes additional 
techniques to optimize both risk and revenue for portfolio 
management. Portfolio optimization (management) normally 
includes two- fold problem, the control of risk which is 

25 usually associated with the volatility (and as it's recently 
understood with higher moments of the multivariate 
distribution of the return of the assets constituting the 
portfolio) . This is a formidable problem and a variety of 
methods have been already proposed to treat it, mention for 

30 instance one -dimensional measures of risk in terms of the 

Value-at-risk (VAR) . Also the problem of the so-called heavy 
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tails has been discussed in the literature in order to manage 
higher order yet less probable risk. 

Mathematically the problem is reduced to 
minimization of the variance and kurtosis of the weighted sum 

5 

y = t V 1 

i = l 



10 of specified stocks under the wealth condition 



i = l 



15 

In what we propose the problem posed is to first 
search for "good" families of stocks on the basis of all 
stocks available at the market. For instance a family of two 
stocks perfectly anticorrelated will simultaneously minimize 
20 variance as well as kurtosis (which defines the higher risk 
contribution) properly created from these stocks. Their 
image in the vector space of zero-mean deviations is 

< > 

25 

and the angle is given by the correlation coefficient 

r(x.x.) 



GOV 



P. 



3 0 XJ ilx.lUlx^l 
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with ||xj being the standard deviation, the norm in that 
space, and a vector 

x ± = jx^ , x? , ... , x* I represents detrended time series of the 
5 stock w i" . 

By calculating the correlation matrix of the stock available 
at the market we can pick up the families consisting of a 
reference vector (stock) w a father" and strongly 

10 anticorrelated with it members "sons" . These members could 
be statistically independent: p ±j = 0 such that they are 
perpendicular to each other. Using standard window 
techniques robustness of the families is checked and most 
persistent members are chosen to constitute the basis from 

15 which optimal portfolios should be built. Minimization of 

the portfolio variance inside the family is achieved by means 
of standard procedure of linear programming or some other 
known technique. As a result, we have a set of "good" 
portfolios suitable to make search in order to optimize and 

20 expected return. 

Given requirements for risk and return the techniques 
described allows to build up corresponding portfolios. The 
whole procedure is comprised of three phases: 

25 • Search for good families. 

a) Creation of the correlation space, including 
shifted correlations . 

b) Looking for xv anticorrelated" families. 

c) Checking robustness of the families. 

30 

• Building perspective portfolios 
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a) Creating two, three and four . stock clusters with 
minimal variance. 

b) Evaluating their risk 

5 • Optimization of expected return (revenue) of a portfolio 
chosen from the perspective ones. 

• The first most difficult' part of the program has been 

checked on the samples consisting of 32 and 100 time 

10 series with intentionally created clusters in the first 

32 of them. There were, four 2 -point clusters, four 3- 
point clusters, and three 4 -point clusters comprising 
all the first 32 series. Applying the method described 
we identified all the clusters in the sample of 32 as 

15 well as in the sample of 100 series. Corresponding 

families are listed in FIG. 27, FIG. 28a and FIG. 28b 
for 

p. . < -0.3 

20 

and on FIG. 29, where are also shown lesser couplings between 
the clusters and some other series. As far as the method is 
geometrical in character, identification of clusters for the 
both samples (32 and 100) only differs in time needed to 
25 build the correlation matrix, it is proportional to n (n-1) /2 . 

Portfolio Risk Management Using Multi-Gaussian Distributions 

Portfolio Risk Management (PRM) is essentially 
3Q based on the techniques applied to estimate probability of 
"bad events" . For that one should know the corresponding 
probability distribution. During several years after Black 
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and Scholes in finance business became widespread the VAR 
machinery which used Gaussian distribution, even though 
market data showed quite noticeable deviation from it. As a 
result, the so-called problem of fat tails arose along with 
5 new fashion to treat it on the ground of Levy distribution 
(later on the truncated Levy distribution) . Underline there 
is no evident reason to regard the latter the only possible 
cure. Moreover even its relation to the problem from the 
market point of view, has still remained foggy. 
10 This invention addresses the problem by using a 

very well-known two-time scale formalism going back to Van 
der Paul and then Bogolubov and Mitropolskii . 

In our case the method gives rise to the cwo- 
gaussian distribution (2-GD) 

15 

x 2 

(1) p» = _^e 2Tl + e 2T > 



2nT. 



2 



20 For the sake of simplicity we consider the symmetric 
distribution with zero mean and the normalization 

(2) «i + «2 = 1 

25 

(3) n^+n^l 

This approach allows us to represent market dynamics at a 
level of description intermediate between macroscopic 
2Q modeling (one-factor models) and microscopic modeling 

(individual agent models or multi- factor models, general 
stochastic processes) . Two-gaussian approach being the 
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simplest: exactly solvable model yet can catch certain 
specific features of the fat tails behavior. It describes 
two groups of investors acting at the market with different 
time scale behavior . 
5 For comparison we are using the standard normalized 

distribution (the Black- Scholes approach) 



x 2 

1 



-e 2 



(4) g( X ) = 

10 



that implies the variance is measured in units of the a. 
Let us consider the probability, sometimes called VAR, of 
15 negative deviations larger than a specified value "a" 

( 5) p (xs-a) =f ~ a p(x)dx 

J -oo 



20 and the difference 



(6) 5P r (x<a) = f 5p(x)dxH J ( Pjfc (x) -g(x) )dx 



25 Since the first two moments are same for the both 

distributions the difference is proportional to the kurtosis 
of 2-GD (1) : 

( 7 ) 5P>.a)^^'3 = 3n A (r 2 - T ± f 

30 
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It can be shown that .under conditions (2) and (3) the 
difference while expanding in the series over moments is 



5(8) 5P =— e 2 — (x 3 -3x) + o(<x 6 >) 



This difference can be incorporated by the so-called implied 
volatility (variance) through considering the effective 
10 (implied) gaussian process 



g(x) - — e 20i 

/2no.. 



im 



15 



where o im = 1 + 5a is the implied variance: 



(9) o. = 1 +— K(a 2 -3) 



20 



In order to check out formula (8) we calculated 5P by using 
the exact distribution (1) and (8) with the kurtosis given by 
(7) . The approximate formula using implied volatility gives 

25 sufficiently good qualitative description and for certain 

values of K and "a" even quantitative one. But there is the 
region of a's where the approximate formula (8) along with 
the notion of implied volatility (9) works bad. Since the 2- 
GD is a very simple model containing enough parameters to fit 

3 q the market data it can be used as an exact tool to calculate 
the "fat tall VAR" . Moreover as far as pricing formulas for 
all the derivatives available at the market linearly depend 
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on the distribution function the model (1) is exact solvable 
in the sense that all results are straightforward. For 
instance the option pricing formula is 



c = e 



" r(r ~ t) f (S T ~ X) g(S T )dS T = 5N(d x ) 'XB' riT t) N(d 2 ) 



where S is the current (stock) price of the underlying asset, 
10 S T is the maturity asset price and X is the strike (exercise) 
price. We see that the price is linear (proportional) over 
the distribution density and in the case of two gaussians it 
will contain four terms with the cumulative gaussian 
functions N(d) instead of two: 



15 



c = e- r(r " t} | (S T - X) [n^lT^S^ + n 2 g [T 2 ; S T ) dS T + n 2 N {T 2 ; d 2 ) ] 

X 

= Sln^T^c^) +n 2 W(T 2 ;d 1 )] -Xe" r(r " fc> [n^iT^d^] 



20 



To make the model less restrictive one can consider three 
(and even more) gaussian distributions 



i = l 

with the normalization conditions 



30 
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E n iTi - 1 

E^J = 1 



such that the difference 

-a. 



can be made proportional to any specified moment. 



Finally we give the next term in the expansion (8) : 




0 with 

<5x 6 > =n 2 (r 2 2 - Tl)[T 2 -l) +K/3 

5 

for the 2-GD. 

FIG. 30 discloses a representative computer system 
3010 in conjunction with which the embodiments of the present 
invention may be implemented. Computer system 3010 may be a 
0 personal computer, workstation, or a larger system such as a 
minicomputer. However, one skilled in the art of computer 
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systems will understand that- the present invention is not 
limited to a particular class or model of computer. 

As shown in FIG. 3 0, representative computer system 
3010 includes a central processing unit (CPU) 3012, a memory 
5 unit 3 014, one or more storage devices 3 016, an input device 
.3 018, an output device 3020, and communication interface 
3022. A system bus 3024 is provided for communications 
between these elements. Computer system 3010 may 
additionally function through use of an operating system such 
10 as Windows, DOS, or UNIX. However, one skilled in the art of 
computer systems will understand that the present invention 
is not limited to a particular configuration or operating 
system. 

Storage devices 3016 may illustratively include one 
15 or more floppy or hard disk drives, CD-ROMs, DVDs, or tapes. 
Input device 3018 comprises a keyboard, mouse, microphone, or 
other similar device. Output device 3020 is a computer 
monitor or any other known computer output device. 
Communication interface 3 022 may be a modem, a network 
20 interface, or other connection to external electronic 
devices, such as a serial or parallel port 

While the above invention has been described with 
reference to certain preferred embodiments, the scope of the 
present invention is not limited to these embodiments. One 
25 skill in the art may find variations of these preferred 

embodiments which, nevertheless, fall within the spirit of 
the present invention, whose scope is defined by the claims 
set forth below. 

30 
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Claims 



1. A method for performing analysis and 
5 prediction of economic markets comprising the steps of: 

initializing a space of economic models ; 

executing one or more of said economic models to 
generate a corresponding one or more predicted time series; 

determining one or more fitnesses corresponding to 
10 said one or more economic models by comparing said plurality 
of predicted time series to actual data; and 

searching for at least one of said economic models 
having a maximal fitness. 

15 2 . A method as in claim 1 further comprising the 

steps of : 

repeating said executing one or more os said 
economic models step, said determining one or more fitnesses 
step, and said searching step to adapt said space of economic 
20 models to changes in the actual data. 

3 . A method for generating models of an economy 
as in claim 1 further comprising the steps of: 

modifying said one or more economic models to 
25 increase said corresponding one or more fitnesses. 

4. A method as in claim 1 further comprising the 

steps of : 

combining at least two of said economic models to 
3 0 create at least one new economic model; and 

adding said at least one new economic model to said 
space of economic models. 
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5. A method as in claim 4 wherein said combining 
step further comprises the step of: 

defining a space of stacking function of said 
economic models; and 
5 applying one or more genetic operators over said 

space of stacking functions to create said at least one 
economic model . 

6. A method as in claim 1 wherein said space of 
10 economic models comprise one or more members of the group 

consisting of AR(n), ARV, ARCH, GARCH, Black Scholes, 
Langevam equations, Levy and truncated Levy. 

7 . A method as in claim 3 wherein said space of 
15 economic models comprises one or more hybrid models wherein 

said hybrid models comprise one or more hybrid models wherein 
said hybrid models comprise one or more components, each of 
said components contributing at least one portion of said 
predicted time series. 

20 

8. A method as in claim 7 wherein said modifying 
said one ore more economic models step comprises the steps 
of: 

adjusting a permutation of said components of said 

25 hybrid models; and 

adjusting said contribution of said components of 

said hybrid models. 

9. A method as in claim 1 wherein said 

30 determining one or more fitnesses step comprises the step of: 
extracting a plurality of observables from said 
predicted time series. 
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10. A method as in claim 9 wherein said 
determining one or more fitnesses step further comprises the 
step of : 

comparing said observables to said actual data. 

5 

11. A method as in claim 9 wherein said extracting 
step uses one or more techniques from the group consisting of 
visual analysis, principal component analysis, data mining, 
localized principal component analysis and independent 

10 component analysis. 

12 . A method as in claim 9 further comprising the 
step of determining a fitness landscape of said one ore more 
fitnesses overs said space of economic models. 

15 

13. A method as in claim 12 wherein said searching 
step comprises the step of : 

searching over said fitness landscape for said at 
least one economic model having the maximal fitness. 

20 

14 . A method as in claim 12 further comprising the 

step of : 

computing at least one probability distribution for 
said plurality of observables. 

25 

15. A method as in claim 14 further comprising the 

step of: 

comparing said probability distribution for said 
plurality of observables to a corresponding probability 
3Q distribution of said actual data; 

determining a distance between one of said economic 
models and another of said economic models ons said fitness 
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landscape by computing a comparison between said probability 
distribution of said one economic model and said probability 
distribution of said another economic model. 

5 16. A method as in claim 15 wherein said 

determining a distance step uses one or more techniques from 
the group consisting of the Randan -Nikodyn Thereom and the 
Kulback-Liebler method. 

10 17. A method as in claim 1 further comprising the 

step of : 

determining at least one state of the economy from 
said actual data. 

15 is. A method as in claim 17 wherein said state of 

the economy comprises a current data and a plurality of past 
data for a plurality of assets xn the economy. 

19. A method as in claim 18 wherein said current 
20 data and said past data comprise prices of said plurality of 

assets . 

20. A method as in claim 17 further comprising the 
step of repeating said determining at least one state of the 

25 economy step over time to assess a measure of stationarity or 
non-stationarity of the economy. 

21. A method as in claim 20 further comprising the 
step of repeating said executing one or more os said economic 

30 models step, said determining one or more fitnesses step, and 
said searching step to adapt said space of economic models to 
changes in said state of the economy. 
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22. A method as in claim. 21 further comprising the 
step of tuning a rate of said searching for at least one of 
said economic models having maximal fitness in accordance 
with said measure of stationarity or non- stationarity of the 

5 economy . 

23. A method as in claim 17 wherein said predicted 
time series are determined from executing said one or more 
economic models in accordance with said state of the economy. 

10 

24. A method as in claim 1 further comprising the 

step of : 

performing risk minimization from said predicted 
time series. 

15 

25. Computer executable software code stored on a 
computer readable medium, the code for performing analysis 
and prediction of economic markets, the code comprising: 

code to initialize a space of economic models ; 
2Q code to execute one or more of said economic* models 

to generate a corresponding one or more predicted time 
series; 

code to determine one or more fitnesses 
corresponding to said one or more economic models by 
25 comparing said plurality of predicted time series to actual 
data; and 

code to search for at least one of said economic 
models having a maximal fitness. 

3Q 26. A programmed computer system for performing 

analysis and prediction of economic markets comprising at 
least one memory having at least one region storing computer 
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executable program code and at least one processor for 
executing the program code stored in said memory, wherein the 
program code comprises : 

code to initialize a space of economic models ; 
5 code to execute one or more of said economic models 

to generate a corresponding one or more predicted time 
series; 

code to determine one or more fitnesses 
corresponding to said one or more economic models by 
10 comparing said plurality of predicted time series to actual 
data ; and 

code to search for at least one of said economic 
models having a maximal fitness. 

15 27. A system for performing analysis and making 

predictions of a securities market comprising 

a model of the securities market comprising: 

a plurality of dealer agents representing a 
corresponding plurality of market makers in the securities 
2 0 market, each of said dealer agents having at least one of a 
plurality of dealer strategies for setting a price of 
securities; 

a plurality of investor agents representing a 
corresponding plurality of investors in the securities 
25 market, wherein said investor agent determines whether to buy 
or sell one or more of the securities at a given price ; 

at least one market model comprising one or 
more rules for exchanging the securities; and 

a simulator for executing said model to perform 
2Q analysis on the securities market. 
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28. A system as in claim 2 7. wherein said plurality 
of dealer agents set a bid price, bid volume, offer price and 
offer volume of a security. 

5 29. A system as in claim 27 wherein said plurality 

of dealer strategies comprise static strategies. 

30. A system as in claim 27 wherein said plurality 
of dealer strategies comprise learning strategies. 

10 

31. A system as in claim 27 wherein said plurality 
of dealer strategies comprise price discovery strategies. 

32. A system as in claim 27 wherein said plurality 
15 of dealer strategies comprise parasitic strategies. 

33. A system as in claim 27 wherein said plurality 
of dealer strategies comprise a basic dealer strategy, said 
basic dealer strategies maintain its price on a security 

2 0 until it receives a predetermined number of trades on the buy 
or sell side wherein it decreases or increases its price on 
the security respectively. 

34. A system as in claim 27 wherein said plurality 
25 of dealer strategies comprise a price volume dealer strategy, 

said price volume dealer strategy deduces demand and supply 
for a security from market data and sets its price for the 
security in accordance with said demand and said supply. 

35. A system as in claim 27 wherein said plurality 
of dealer strategies comprises a new volume dealer strategy 
wherein said new volume dealer strategy, said new volume 
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dealer strategy raises or lowers its price on the security 
when it observes more buys or more sells respectively for the 
security. 

5 36. A system as in claim 27 wherein said plurality 

of dealer strategies comprise a parasitic dealer strategy, 
said parasitic dealer strategy waits for a sufficiently 
narrow spread with sufficient volume for a security before 
buying or selling the security inside the spread. 

10 

37. A system as in claim 27 wherein said plurality 
of dealer strategies comprises a new parasite strategy, said 
new parasite strategy decides whether to buy or sell the 
security inside the spread based on observed excess supply 
15 and/or excess demand for the security. 



38. A system as in claim 27 wherein said plurality 
of dealer strategies comprise a matching dealer strategy, 
said matching dealer strategy learns connections between one 
2 0 ore more observations and one or more actions and makes 
decisions based on said learned connections. 



39. A system as in 
observations comprises one or 
25 inventory level, spread size, 
value and trade volume. 



claim 38 wherein said 

more of the group consisting of 

bid price relative to market 



40. A system as in claim 38 wherein said actions 
comprises one or more of the group consisting of raising a 
3q price, increasing a bid and narrowing the spread. 
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41. A system as in claim 27 wherein said plurality 
of dealer strategies comprise a classifier dealer strategy, 
said classifier dealer strategy learns pattern connections 
between one or more patterns of observations and one or more 

5 actions and makes decisions based on said learned pattern 
connections . 

42. A system as in claim 2 7 wherein said plurality 
of dealer strategies comprise a dynamical system dealer, said 

10 dynamical system dealer uses a discrete dynamical system to 
set its bid price, ask price and mean price. 

43. A system as in claim 27 wherein said plurality 
of investor agents comprises a degree of informedness 

15 representing a measure of said investor agents knowledge of a 
true value of a security. 

44 . A system as in claim 27 wherein said one or 
more rule for exchanging securities comprises at least one 

20 tick price . 

45. . A system as in claim 27 further comprising a 
plurality of parameters that quantify an ability of the 
market to track prices of the securities. 

25 

46. A system as in claim 45 wherein said 
parameters comprises one or more members of the group 
consisting of the root mean square of the difference between 
a bid price and a time value, a root mean square of the 

30 difference between an ask price and a true value, an average 
value of a spread and a standard deviation of the spread. 
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47. Computer executable software code stored on a 
computer readable medium, the code for performing analysis 
and making predictions of a securities market, the code 
comprising: 

5 code to model the securities market comprising: 

code to represent a plurality of market makers 
in the securities market including a plurality of dealer 
strategies for setting a price of securities; 

code to represent a plurality of investors in 
10 the securities market and to determine whether to buy or sell 
one or more of the securities at a given price ; 

code to represent the market including one or 
more rules for exchanging the securities; and 

code to simulate the model to perform analysis on 
15 the securities market . 

48. A programmed computer system for performing 
analysis and making predictions of a securities market 
comprising at least one memory having at least one region 

2 0 storing computer executable program code and at least one 
processor for executing the program code stored in said 
memory, wherein the program code comprises: 

code to model the securities market comprising: 

code to represent a plurality of market makers 
25 in the securities market including a plurality of dealer 
strategies for setting a price of securities; 

code to represent a plurality of investors in 
the securities market and to determine whether to buy or sell 
one or more of the securities at a given price; 
3Q code to represent the market including one or 

more rules for exchanging the securities; and 
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code to simulate the model to perform analysis on 
the securities market. 

49. A method for creating a portfolio of assets 
5 comprising the steps of: 

initializing at lease one current portfolio of 
assets having one or more desired properties; 

generating at least one next portfolio of assets to 
form said current portfolio, said next portfolio having a 
10 smaller size then said current portfolio; 

repeating said generating step after replacing said 
current portfolio with said next portfolio while said next 
portfolio has the one or more desired properties; 

defining a size range having an upper limit equal 
15 to the size of said current portfolio and a lower limit equal 
to the size of said next portfolio; and 

searching for an optimal portfolio of assets to 
within said size range, 

20 50. A method as in claim 49 wherein said 

generating at least one next portfolio comprises the steps 
of: 

sampling said assets from said at least one current 
portfolio to form a plurality of candidate portfolios; and 
25 selecting said at least one portfolio from said 

candidate portfolios . 

51. A method as in claim 50 wherein said selecting 
said at least one next portfolio step comprises the step of 
30 choosing said candidate portfolios that are optimal with 
respect to the one or more desired properties. 
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52. A method as in claim 50. wherein said selecting 
said at least one next portfolio step comprises the step of: 

choosing at least one of said candidate portfolios 
that is pareto optimal with respect to the one or more 
5 desired properties. 

53. A method as in claim 50 wherein said sampling 
of said assets has a uniform distribution. 

10 54. A method as in claim 50 wherein said sampling 

of said assets has a distribution that varies with time. 

55. A method as in claim 54 wherein said sampling 
step further comprises the steps of : 

15 altering said distribution over time based on said 

one or more desired properties over previous choices of said 
candidate portfolios. 

56. A method as in claim 55 wherein said altering 
20 said distribution over time step boosting. 

57. A method as in claim 49 wherein said searching 
for an optimal portfolio of assets step comprises the steps 
of: 

25 determining an intermediate size within said size 

range ; 

generating a lease on intermediate portfolio of 
assets of said intermediate size from said current portfolio; 
and 

3Q if said intermediate portfolio has the one or more 

discreet properties 
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modifying said upper limit of said size range 
to equal said intermediate size; 
else 

modifying said lower limit of said size range 
5 to equal said intermediate size. 

58. A method as in claim 57 wherein said searching 
for an optional portfolio of assets to step further comprises 
the step of: 

10 repeating said steps of claim 53 to converge to 

said optimal portfolio of assets. 

59. A method as in claim 49 wherein said size of 
said current portfolio is defined as the number of said 

15 assets in said current portfolio. 

60. A method as in claim 49 wherein said size of 
said next portfolio is a fraction of said size of said 
current portfolio. 



20 



61 . 



A method as in claim 60 wherein said fraction 



is M- 



62 . 



A method as in claim 60 wherein said 



25 generating at least one next portfolio of assets step 
comprises the step of: 

tuning said size of said next portfolio. 



63 . 



A method as in claim 4 9 further comprising the 



3Q step of: 



determining the number of said at least one next 



portfolios to generate. 
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64 . 



A method that as in claim 63 form the 



comprising the step of timing said number of said at least on 
next portfolios to generate. 



5 



65 . 



A method as in claim 4 9 wherein said one or 



more desired properties comprise a low value at risk. 

66. Computer executable software code stored on a 
computer readable medium, the code for creating a portfolio 
10 of assets, the code comprising: 

code to initialize at lease one current portfolio- 
of assets having one or more desired properties; 

code to generate at least one next portfolio of 
assets to form said current portfolio, said next portfolio 
15 having a smaller size then said current portfolio; 

code to repeat said generating step after replacing 
said current portfolio with said next portfolio while said 
next portfolio has the one or more desired properties; 

code to define a size range having an upper limit 
20 equal to the size of said current portfolio and a lower limit 
equal to the size of said next portfolio; and 

code to search for an optimal portfolio of assets 
to within said size range. 



portfolio of assets comprising at least one memory having at 
least one region storing computer executable program code and 
at least one processor for executing the program code stored 
in said memory, wherein the program code comprises: 

code to initialize at lease one current portfolio 
of assets having one or more desired properties; 



25 



67 . 



A programmed computer system for creating a 
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code to generate at least one next portfolio of 
assets to form said current portfolio, said next portfolio 
having a smaller size then said current portfolio; 

code to repeat said generating step after replacing 
5 said current portfolio with said next portfolio while said 
next portfolio has the one or more desired properties; 

code no define a size range having an upper limit 
equal to the size of said current portfolio and a lower limit 
equal to the size of said next portfolio; and 
10 code to search for an optimal portfolio of assets 

to within said size range. 

68. A method for optimizing a portfolio of assets 
comprising the steps of ; 

15 determining a fitness landscape representation over 

a space of portfolios; 

determining at least one optimal searching distance 
in said fitness landscape representation; and 

searching for optional ones of said portfolio of 
20 assets at said optimal distances. 

69. A method as in claim 68 wherein said portfolio 
of assets is defined as a vector corresponding to said assets 
wherein each element of said vector identifies the umber of 

25 said corresponding assets in said portfolio. 

70. A method as in claim 69 wherein said distances 
between a first of said portfolio and the second one of said 
portfolios is defined as the difference between said vector 

2q of said first portfolio and said vector of said second 
portfolio. 
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71. A method as in claim 68. wherein said fitness 
of said landscape representation comprises a value at risk. 

72 . A method as in claim 68 wherein said 

5 determining a fitness landscape representation step comprises 
the step of inferring said fitness landscape from historical 
data. 

73 . Computer executable software code stored on a 
10 computer readable medium, the code for optimizing a portfolio 

of assets, the code comprising: 

code to determine a fitness landscape 
representation over a space of portfolios; 

code to determine at least, one optimal searching 
15 distance in said fitness landscape representation; and 

code to search for optional ones of said portfolio 
of assets at said optimal distances. 

74 . A programmed computer system for optimising a 
20 portfolio of assets comprising at least one memory having at 

least one region storing computer executable program code and 
at least one processor for executing the program code stored 
in said memory, wherein the program code comprises: 

code to determine a fitness landscape 
25 representation over a space of portfolios; 

code to determine at least one optimal searching 
distance in said fitness landscape representation; and 

code to search for optional ones of said portfolio 
of assets at said optimal distances. 

30 

75. A method for creating a portfolio of a 
plurality of assets comprising the steps of: 
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determining a plurality of anti -correlated families 
wherein each of said ant i -correlated families contains two or 
more of said plurality of assets that are anti -correlated ; 
and 

5 determining one or a more perspective portfolios of 

said plurality of assets from said plurality of anti- 
correlated families . 



76. A method in claim 75 further comprising the 

10 step of : 

receiving one or more time series corresponding to 
one or more of said plurality of assets. 

77. A method as in claim 76 further comprising the 

15 step of : 

determining correlations among said plurality of 

assets . 

78. A method as in claim 75 wherein said 

20 determining a plurality of anti -correlated families step 
comprises the steps of: 

selecting at least one of said plurality of assets 
to be a first member of one of said ant i -correlated families; 

determining one or more of said plurality of assets 
25 having a value of said correlation with said first member of 
said ant i -correlated family that is less than a negative 
limit; and 

including said one or more assets having a 
correlation value that is less than the limit as additional 
3q members in said anti -correlated family. 
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79. A method as in claim 78 further comprising the 
step of separating said steps of claim 4 for another one of 
said plurality of assets to determine another of said 
plurality of ant i -correlated families. 

5 

80. A method as in claim 77 further comprising the 

step of: 

eliminating drift from said one or more time series 
to create corresponding one or more detrended time series. 

10 

81. A method as in claim 80 wherein said 
correlations among said plurality of assets is determined 
from said one or more time series. 

15 82. A method as in claim 77 wherein said 

correlations are shifted in time. 

83. A method as in claim 75 wherein said 
determining one or more perspective portfolios step comprises 

20 the step of selecting one of said plurality of anti- 
correlated families ; 

determining percentages of said assets from said 
selected family to include in one of said perspective 
portfolios. 

25 

84. A method as in claim 83 further comprising the 
steps of repeating said selecting step and said determining a 
percentage of each of said assets step with different ones of 
said anti-correlated families to determine said one or more 

3q perspective portfolios. 
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85. A method as in claim 83 wherein said 
percentages of said assets from said selected family are 
determined from said perspective portfolio with a minimal 
variance . 

5 

86. A method as in claim 75 further comprising the 
step of choosing at least one of said perspective portfolios 
as said portfolios of assets. 

10 87 . A method as in claim 86 for the comprising the 

step of optimizing an expected return of said portfolio of 
assets . 

88. Computer executable software code stored on a 
15 computer readable medium, the code for creating a portfolio 
of a plurality of assets, the code comprising: 

code to determine a plurality of anti -correlated 
families wherein each of said anti -correlated families 
contains two or more of said plurality of assets that are 
20 anti-correlated; and 

code to determine one or a more perspective 
portfolios of said plurality of assets from said plurality of 
anti-correlated families. 



25 8 9. A programmed computer system for creating a 

portfolio of a plurality of assets comprising at least one 
memory having at least one region storing computer executable 
program code and at least one processor for executing the 
program code stored in said memory, wherein the program code 

30 comprises: 

code to determine a plurality of anti-correlated 
families wherein each of said anti-correlated families 
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contains two or more of said plurality of assets that are 
ant i -correlated ; and 

code to determine one or a more perspective 
portfolios of said plurality of assets from said plurality of 
5 anti -correlated families. 

90. A method of modeling a portfolio of assets 
comprising the step of; 

defining a probability distribution of the value of 
10 the portfolio of assets with two or more Gaussian 
distributions. 

91.. A method as in claim 90 wherein said two or 
more Gaussian distributions comprise at least one parameter. 

15 

92. A method as in claim 91 further comprising the 

step of : 

receiving historical data on the value of the 
portfolio of assets; and 
20 determining values for said at least one parameter 

from said historical data. 

93 . A method as in claim 91 wherein said defining 
a probability distribution a probability distribution step 

25 comprises the step of: 

defining said probability distribution as a sum of 
said two or more Gaussian distributions. 

94 . A method as in claim 93 wherein said defining 
30 said probability distribution as a sum step comprises the 

step of defining said probability distribution as: 
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«7«< n i' r i> = £n.g(l,r.) 

n t represent normalisation parameters 
T represent variance ■ 

5 

95. A method that as in claim 94 wherein said 
defining said probability distribution as a sum step further 
comprises the step of: 

10 defining a plurality of normalization conditions 

as : 

En, - 1 
Zn~ Ti = 1 

En, 7f = 1 

15 

96. A method as in claim 90 for further comprising 
the step of : 

computing a value at risk from said probability 
distribution of said two or more Gaussian distributions . 

20 

97 . A method as in claim 90 for further comprising 
the step of : 

computing an option price from said probability 
distribution of said two or more Gaussian distributions. 

25 

98. A method as in claim 97 wherein said computing 
an options price step comprises the steps of: 

determining an option pricing formula as: 

c = e- r(T - fc) f (S T - X)g(S t )dS T 

30 

= S[n 1 N(T 1 ;d 1 ) + n 2 N ( T 2 ; d x ) } - Xe ~ rlT ' t} [n^iT^d^] 
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wherein : 

T represents variance 

N represents normal density function 
5 t represents time 

99. Computer executable software code stored on a 
computer readable medium, the code for modeling a portfolio 
of assets , the code comprising : 

10 code to define a probability distribution of the 

value of the portfolio of assets with two or more Gaussian 
distributions . 

100. A programmed computer system for modeling a 
15 portfolio of assets comprising at least one memory having at 

l eas t one region storing computer executable program code and 
at least one processor for executing the program code stored 
in said memory, wherein the program code comprises: 

code to define a probability distribution of the 
20 value of the portfolio of assets with two or more Gaussian 
distributions . 



25 
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D. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Chicago 


STATE OR FOREIGN COUNTRY 

IL 


COUNTRY OF CmZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

5448 Ridgewood Court 


cm- 
Chicago 


STATE OR COUNTRY 
IL 


ZIP CODE 

60615 






SIGNATURE OF INVENTOR 204 ^ 

3«s e/e P. Cavo^^ 


•% 

DATE* - 


2 
0 
5 


FULL NAME 
OF INVENTOR 


LAST NAME 

TEILTELBAUM 


FIRST NAME 

Daniel 


MIDDLE NAME 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 

NM 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

618 Old Santa Fe Trail #G 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 205 


DATE 



I hereby declare -that all statements made herein of my own knowledge are true and that all statements made on information and belief are 
believed to be true; and further that these statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements 
may jeopardize the validity of the application or any patent issuing thereon. 




FULL NAME 
OF INVENTOR 


LAST NAME 

OUTKIN 


FIRST NAME 

Alexander 


MIDDLE NAME 


2 
0 
6 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

Russia 




POST OFFICE 
ADDRESS 


STREET 

1548 Avida de Las Americas 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 20 1 


DATE 

3/2<?/l?{ 




FULL NAME 
OF INVENTOR 


LAST NAME 

SAWHILL 


FIRST NAME 

Bruce 


MIDDLE NAME 

K. 


2 
0 
7 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

United States 




POST OFFICE 
ADDRESS 


STREET 

1028 Don Cubero 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 202 / n 




DATE 










f / zs/z&o f 




FULL NAME 
OF INVENTOR 


LAST NAME 

MAKHANKOV 


FIRST NAME 

Vladimir 


MIDDLE NAME 

G. 


2 
0 
8 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

Russia 




POST OFFICE 
ADDRESS 


STREET 

4148 Chaparron Place 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 203 


DATE 






FULL NAME 
OF INVENTOR 


LAST NAME 

DARLEY 


FIRST NAME 

Vince 


MIDDLE NAME 
M. 


2 
0 
9 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

United Kingdom 




POST OFFICE 
ADDRESS 


STREET 

794 Calle Altamira 


CITY 

Santa Fe 


STATE OR COUNTRY 

Nrv£ 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 204 




DATE 
















FULL NAME 
OF INVENTOR 


LAST NAME 


FIRST NAME 


MIDDLE NAME 


2 
0 
10 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 




POST OFFICE 
ADDRESS 


STREET 


CITY 


STATE OR COUNTRY 


ZIP CODE 




SIGNATURE OF INVENTOR 205 


DATE 





. CJ-J 



I hexxby declare th« ail sttccmcriis made herein of my own JcnowJcdgc are true and thai all rUCcmcou made on information and belief arc 
bcHcvcd to be true; And further thai these statements were made wftb Ac knowledge thai willful Gdsc jtucmcnte arid vho like so made arc 
punishable by flnc or imprisonment, or botfi, under Section 1 00 1 of Title 18 of ine United Statas Code and that cuch willful felsc statements 
mayJ«opardix« the vaJfdlty of the application or any paUnt tosulnt thac«oa, 



2 
0 
1 


FULL NAME 
OF INVENTOR 


KAUFFMAN 


TOUT NAME 

Stuart 


MIDOUMAUS 

A. 


RESIDENCE A 
CniZENSHIP 


COY 
S&MH Fc 


iunon wnsicu country 
NM 


COUNTRY Of CTTCKUl KJf 

United Suites 


POST OFFICE 
ADDRESS 


misr 

1811 S. Curtid Cruz Blanca 


/ Santa Fa 


STATE 0* COCNTKY 

NM 


xirccoe 
67501 






OATI 

<={*Z f'ZG 




2 
0 
2 


FUULKAME 
OF INVENTOR 


LAJTNAMH ^ 

. MAC READY 


William 


a 


RESIDENCE & 
CITIZENSHIP 


crry 

_£«ita-Ea — 


fTATK OA POIUOCW COtJ»m.Y 


cownv or CmZXNCMlr 

Canada 


POST OFFICE 
ADDRESS 




cmr 
Santa Fc 


NM 


CP CODC 




liOMATunf or rv/vorro*. 303 


OaTH 

tec tr 




2 
0 
3 


FULL NAME 
OF INVENTOR 


PALMER 


raj r wax* 
Richard 


G. 


RESIDENCE 4 
CITIZENSHIP 


cmr 

Durham 


rTATl OR FDBJ1CM CfttVTav 

NC 


cowir-r 0 f arxuhSHT 
Untied Kingdom 


POST OFFICE 
ADDRESS 


4030 Forcstdalc Drive 


CTTY 

Durham 


STATCOR.COUKTHY 

NC 


flf coot 
27712 




sk^atuh£ or wvo/roa it»i 


tMTtt 




FULL NAME 
OF INVENTOR 


COWAN 


PUlTNAMf 

Jai* 


MIOOLSVAAa 

D. 


0 
4 


RESIDENCE 8. 

crrizENSHtP 


cmr 
CMcqgo 


STATE OR P0*£lCw COUNTRY 
IL 


United States 




POST OFFICE 
ADDRESS . 


344$ Ridgewood Court 


cmr 
Chicago 


rTATtoacowntY 
IL 


ZftCODfi 

$0615 




OF IWVSMT013U4 


* * 




2 " 

0 
5 


FULL NAME 
OF INVENTOR 




Daniel 


MUWKMAMr 


RESIDENCE A 
CmZENSHZP 


CTTY 


STATU 0€ M>«j OX OQITHTT / A 


count* v or cmeswHir 
United Sww 




POST OFFICE 
ADDRESS 


OM Guuu Pl T)tfl #0 


arr ty^cr^ 
.8*161 To 




tircooe 




jUtt^^f * 


DATS 





I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are 
believed to be true;* and further that these statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements 
may jeopardize the validity of the application or any patent issuing thereon. 





FULL NAME 
OF INVENTOR 


LASTNAME 
1/ AT IPPN/t A M 


FIRST NAME 
C( llnr ( 


MIDDLE NAME 

A. 


2 
0 
1 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
fNlVl 


COUNTRY OF CITIZENSHIP 

United States 




POST OFFICE 
ADDRESS 


STREET 

1811 S. Ojm^no Cruz Blanca 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATURE CffJNVENTOR 201 ^/^Z- 


^ 


DATE 










<f- — - 






FULL NAME 
OF INVENTOR 


LAST NAME 

MACREADY 


FIRST NAME 

William 


MIDDLE NAME 

G. 


2 
0 
2 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

Canada 




POST OFFICE 
ADDRESS 


STREET 

339 1 / 3 -Delgado 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 202 


DATE 




FULL NAME 
OF INVENTOR 


LAST NAME 

PALMER 


FIRST NAME 

Richard 


MIDDLE NAME 

G. 


2 
0 
3 


RESIDENCE & 
CITIZENSHIP 


CITY 

Durham 


STATE OR FOREIGN COUNTRY 
NC 


COUNTRY OF CITIZENSHIP 

United Kingdom 




POST OFFICE 
ADDRESS 


STREET 

4030 Forestdale Drive 


CITY 

Durham 


STATE OR COUNTRY 

NC 


ZIP CODE 

27712 






SIGNATURE OF INVENTOR 203 


DATE 




FULL NAME 
OF TNVFNTDR 


LAST NAME 

COWAN 


FIRST NAME 

Jack 


MIDDLE NAME 

D. 


2 

0 
4 


RESIDENCE & 
CITI7FNSHIP 


CITY 

Chicago 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 

United States 




POST OFFICE ~~ 


—STREET ' " " 

5448 Ridgewood Court 


CITY 

Chicago 


STATE OR COUNTRY 

IL 


ZIP CODE 

60615 






SIGNATURE OF INVENTOR 204 ^ 

^ 3« elc ^. G»*j-j^ 


-» 

DATE 4 - * 




FULL NAME 
OF INVENTOR 


LAST NAME 

TEILTELBAUM 


FIRST NAME 

Daniel 


MIDDLE NAME 


2 
0 
5 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

United States 




POST OFFICE 
ADDRESS 


STREET 

618 Old Santa Fe Trail ftG 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 205 


DATE 



(2) 



DCl - 300946.1 



IN THE UNITED STATES FATRNT AND TRADEMARK OFFICE 
Application of: Stuart A. KAUFFMAN cl ul. 



U.S. Applicaiion No.: 09/868,9^1 Group Art Unit: To be assigned 

U.5L Application Viling Date', June 22, 2001 

Tntaninlional Application Mo.: PCT/US 99/3 0640 Examiner: To be assigned 
International Application Filing Dale: 
December 22, 1999 

For: A S YSTF M AND MKTHOD FOR THE Attorney Docket No.: 9392-0 1 6-990 
ANALYSIS AND PREDICTION OF 
ECONOMIC MA1UCBTS 

CONSENT BY TH E AS&TGNKF UNDER 37 CF,R. 6 Mffldj 

Assistant Commissioner for Patents 
Washington* D.C. 20231 

Sir; 

Bios Group, Inc., having a place of business at 317 Pasco dc Perdta, Santa 
Fe, Now Mexico 87501 , as Assignee of the entire right, title and interest in, to and under the 
U.S. application filed on Juno 22, 2001 entitled "A SYSTEM AND MBTHOD FOR THE 
ANALYSIS AND PREDICTION OF ECONOMIC MARKETS" hereby consents to the 
addition of Vince M. Darlcy &s a named inventor in US. application no. 09/868,98 1, 

Bios Group > Inc., certifies that Jt Is the Assignee of the entire right, title and 
Interest in, to and under U,S. Application, 09/868,981 pursuant to the assignment executed 
by Stuart A. Kauftmon, William G, Macready, Richard G. Palmer, Jack D. Cowan, Daniel 
Teitclbaum, Alexander Outkin, Bruce K. Sawhill, Vladimir G, Makhonkov and Vinco M. 
Darlcy. Documentary evidence of a chain of title trom the original owner to the assignee is 
indicated by the enclosed copy of the assignment. 

I hereby declare thet all statements made herein of my own knowledge aire 
true and that ail statements made on information and belief are believed to be true; and 

DC! -1)2082 i 



further ihiU these statements were made with the knowledge that willful false slateraenis and 
the like so mad* are punishable by fine or unprLionmcnt, or both, under section 1001 of 
Title 18 of the United States Code and that such willful false statements may jeopardize the 
validity of the Patent. 



Signed this day of ^ 



2001. 



Bios Group, Inc. 




Robert R. MacDonald 
Title; President 



DC1 



JN THE UNITED STATES PATENT AND TRADEMARK OFFICE 
Application of; Stuart A, KAUFFMAN et al 

U.S. Application No. : 09/868,98 \ Group Ait Unit: To be assigned 

U.S. Application Filing Date; Juno 22, 2001 

lntcrnalioniil Application No.: PCT/US99/30640 Examiner: Ta be assigned 
International Application Filing Dale: 
Decembfcr 22, 1999 

For: A SYSTEM AND METHOD FOR THE Attorney Docket No.; 9392-016-999 
ANALYSIS AND PREDICTION OF 
ECONOMIC MARKETS 

STATEMENT OF VfNCE M. PARLEY UNDER 37 CF.H 6 t .497^ 

Assistant Commissioner for Patents 
Washington, D.C 20231 

Sir; 

1, Vinco M. Darloy, hereby state that I was mistakenly not included sus a 
named inventor of international application US99/30640 and that this error in Inventorship 
in this international application occurred without any dtsccptive intention on iny part. 

I hereby declare that all statements madfi herein o ( my own knowledge are 
true and that all statements made on information imd belief arc believed to be true, and 
further that these statements were nude with the Icnowledge that willful false statements and 
the like so made arc punistablo by fino or Imprisonment or both, under section 1001 of Tizlc 
18 of the United S tines Code and that such wilifu] false statement may jeopardize the 
validity of (he patent 

Vincc M. Darley 




ASSIGNMENT 



WHEREAS, WE, 

Stuart A. KAUFFMAN, a citizen of the United States, residing at 1811 S. Cam ino Cruz Blanca Santa Fe 
New Mexico 87501; 

William G. MACREADY, a citizen of Canada, residing at 3391/2 Delgado, Santa Fe, New Mexico 87501; 
Richard G. PALMER, a citizen of the United Kingdom, residing at 4030 Forestdale Drive, Durham, North CarolL 



Jack D. COWAN, a citizen of the United States, residing at 5448 Ridgewood Court, Chicago, Illinois 606 15; 

Daniel TEI£TELBAUM, a citizen of the United States, residing at 618 Old Santa Fe Trail UG Santa Fe New 
Mexico 87501; fat£l$W*\ ' 

Alexander OUTKIN, a citizen of Russia, residing at 1548 Avida de Las Americas, Santa Fe, New Mexico 8750 1 ; 

Bruce K. SAWHILL, a citizen of the United States, residing at 1028 Don Cubero, Santa Fe, New Mexico 



Vladimir G. MAKHANKO V, a citizen of Russia , residing at 4 1 48 Chaparron Place, Santa Fe New Mexico 8750 1 - 
and ' 

Vince M. DARLEY, a citizen of the United Kingdom, residing at 794 Calie Altamira, Santa Fe New 
Mexico 87501; 



are the inventors of the invention in A SYSTEM AND METHOD FOR THE ANALYSIS AND PREDICTION OF 
ECONOMIC MARKETS for which we have executed an application for a Patent of the United States 

□ which is executed on □ even date herewith or n _____ 



H which is identified by Pennie & Edmonds llp docket no. 9392-0 16-999 ; 
h which was filed on June 22, 200 1 , Application No. 09/868,98 1 

and WHEREAS, Bios Group, Inc., having a business address of 317 Paseo de Peralta, Santa Fe, New Mexico 87501 
ASSIGNEE is desirous of obtaining our entire right, title and interest in, to and under the said invention and the said 
application: 

NOW, THEREFORE, in consideration of the sum of One Dollar ($1.00) to us in hand paid, and other good and 
valuable consideration, the receipt of which is hereby acknowledged, we, the said ASSIGNORS, have sold, assigned 
transferred and set over, and by these presents do hereby sell, assign, transfer and set over, unto the said ASSIGNEE its ' 
successors, legal representatives and assigns, our entire right, title and interest in, to and under the said invention, and the said 
United States application and all divisions, renewals and continuations thereof, and all Patents of the United States which may 
be granted thereon and all reissues and extensions thereof; and all applications for industrial property protection, including, 
without limitation, all applications for patents, utility models, and designs which may hereafter be filed for said invention in 
any country or countries foreign to the United States, together with the right to file such applications and the right to claim for 
the same the priority rights derived from said United States application under the Patent Laws of the United States, the 
International Convention for the Protection of Industrial Property, or any other international agreement or the domestic laws 
of the country in which any such application is filed, as may be applicable; and all forms of industrial property protection, 
including, without limitation, patents, utility models, inventors' certificates and designs which may be granted for said invention 
in any country or countries foreign to the United States and all extensions, renewals and reissues thereof; 

AND WE HEREBY authorize and request the Commissioner of Patents and Trademarks of the United States, and 
any Official of any country or countries foreign to the United States, whose duty it is to issue patents or other evidence or forms 
of industrial property protection on applications as aforesaid, to issue the same to the said ASSIGNEE, its successors, legal 
representatives and assigns, in accordance with the terms of this instrument. 

AND WE HEREBY covenant and agree that we have full right to convey the entire interest herein assigned, and that 
we have not executed, and will not execute, any agreement in conflict herewith. 

AND WE HEREBY further covenant and agree that we will communicate to the said ASSIGNEE, its successors, 
legal representatives and assigns, any facts known to us respecting said invention, and testify in any legal proceeding, sign all 
lawful papers, execute all divisional, continuing, reissue and foreign applications, make all rightful oaths, and generally do 
everything possible to aid the said ASSIGNEE, its successors, legal representatives and assigns, to obtain and enforce proper 
protection for said invention in all countries. 



signatures. 



Date G - Z^F ~ c? / 



State of ) 
) SS.: 

County of ) 

In the State of /M fr,^ t j-igv^co county of n ^ rvy , on ±1 &fT z<j&[ , before 

me, Srg'jg^, Notary Public, personally appeared Stuart A. KAUFFMANL personally known to me on the basis 

of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 




State of ) 
) SS.: 

County of ) 

In the State of 6 county of "iSwaj rW , on // - /£ -o/ } before 

me, l)fVL^v5 '$7ttt£- , Notary Public, personally appeared William G. MACREADY, personally known to me on the basis 
of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 




WITNESS my hand and official seal 




Date 



f 0/g/c7/ 



2001 



L.S. 



Richard G.. PALMER 



State of N*C 



) SS.: 



County of T>ufttWM ) 

county of 



In the State of < 



on 



(Qd.J^ol, before 



me IT>6A/A//? LtRgaeQ.. Notary Public, personally appeared Richard G. PALMER, personally known to me on the basis 
of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 




Mw frumnminfl ti$x*rfflJ. £0.Z6Qt£ 



Datej . , 2l L.S. 

Jack D. COWAN 

State of ) 
) SS.: 

County of ) 

In the State of , county of on , before 

me, , Notary Public, personally appeared Jack D. COWAN, personally known to me on the basis of 

satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



Date , 200 1 [ L S 

Daniel TEILTELBAUM 

State of ) 
) SS.: 

County of ) 

In the State of , county of , on , before 

me, ^ Notary Public, personally appeared Daniel TEILTELBAUM. personally known to me on the basis 

of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



Date Jfy^ frsfaU r 

7 [ Alexander OUTKIN 

State of ) 
) SS.: 

County of ) 

In the State of H-Xei^ k^co , county of 5^°^ ^fc , on 9- 2 S - <^>{ , before 

me /DA y^&tP . Notary Public, personally appeared Alexander OUTKIN. personally known to me on the basis of 
satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 




Datc_ 




State of ) 
) SS.: 

County of ) 

In the State of f^e^> ^e*o <-<=> , county of ^A^fA , on j3 "Z&QK before 

me. Twi.a5 SVee^g" Notary Public, personally appeared Bruce K. SAWHILL , personally known to me on the basis 
of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 




L.S. 



Vladimir G. MAKHANKOV 



State of / ) 

) SS.: 

County of ) 

In the State of ^e*^ **e.><u county of 1>*^tVw V<^- , on //- , before 

me/ PAi^A*; 6-rfce<^- , Notary Public, personally appeared Vladimir G. MAKHANKOV, personally known to me on the 
basis of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me 
that he executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon 
behalf of which the person acted, executed the instrument. 



WITNESS my hand and official seal 




Date , 200 1 _ L.S. 

Vince M DARLEY 

State of ) 
) SS.: 

County of ) 

In the State of , county of on , before 

me, , Notary Public, personally appeared Vince M. DARLEY, personally known to me on the basis of 

satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. v 



WITNESS my hand and official seal 



State of 
County of 

In the State of 



) SS.: 
) 



3 county of *6ftA^-n\ ^ g 



, on IS **?7 before 



me Tw^m e ^ Notary Public, personally appeared Bruce K. SAWHILL, personally known to me on the basis 
of saSactorv evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 

WITNESS my hand and official seal 




Date 



_,2001 



L.S. 



Vladimir G. MAKHANKOV 



State of 
County of 



In the State of. 



) SS.: 
) 



county of 



, on 



, before 



, ^UUlll) Ul 5 — 7 

me Notary Public personally appeared Vladimir G. MAKHANKOV, personally known to me on the 

bas'is of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me 
that he executed the same in his authorized capacity, and that by his signature on the mstrument the person, or the entity upon 
behalf of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



Date 



JO 



2001 




Vince M. DARLEY 



L.S. 



State of ) 

) SS.: c 

County of ) « 

In the State of , county of , on before 

— Notary Public, personally appeared Vince M. DARLEY, personally known to me on the basis of 
satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed I {he same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed die instrument. 

WITNESS my hand and official seal 



< 



,o o ;;!;:;?• .j:;]; ;o -o 



Date 



State of ) 
) SS.: 

County of ) 



_____ L.S. 

Jack D. COWAN 



In the State of , county of , on , before 

mCj Notary Public, personally appeared Jack D. COWAN, personally known to me on the basis ot 

satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



Dat e iV'Qj'l 



,,2001 



Hj to- 



L,S. 



State of 
County of 



Daniel TEI^TELB AUM 



) SS.: 
) 



re 



/o * *t * o { , before 



In the State of ^&vo c o T county of , 

mc ^A^Ar^rg^c^ Notary Public, personally appeared Daniel TEIKTELB AUM, personally known to me on the basis 
of satisfactory evidence to be the person whose name is subscribed to the" within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument 



WITNESS my hand and official seal 




L.S. 



State of 
County of 



) SS.: 
) 



In the State of. 



Alexander OUTKIN 



county of _ 



, before 



me, , Notary Public, personally appeared Alexander OUTKIN, personally known to me on the basis of 

satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



/ 



Date 



Jack D. COWAN 



State of /U-Z/dO/S 




L.S. 



County of (LOOK- ) 



) SS.: 



In the State of 



_, county of doo/ZL. 



, on 



, . , // -bo -Of , before 

me, ^ATl^tCJ^dJH^^C^UjC f^oXsxs Public, personally appeared Jack D. C OWAN, personally known to me on the basis of 
satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 





•OFFICIAL SEAL* 
Patricia K Chancellor 
Kotaiy Public, State of IHinofs 
My Commission Expires 9-25-2002 



Date , 200 1 ' L.S. 

Daniel TEILTELBAUM 

State of ) 
) SS.: 

County of ) 

In the State of , county of , on , before 

me, . , Notary Public, personally appeared Daniel TEILTELBAUM, personally known to me on the basis 

of satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



Date * /lS /*! , 2001 ' L.S. 

7 Alexander OUTKIN 

State of ) 
) SS.: 

County of ) 

In the State of , county of , on , before 

me, , Notary Public, personally appeared Alexander OUTKIN, personally known to me on the basis of 

satisfactory evidence to be the person whose name is subscribed to the within instrument and acknowledged to me that he 
executed the same in his authorized capacity, and that by his signature on the instrument the person, or the entity upon behalf 
of which the person acted, executed the instrument. 



WITNESS my hand and official seal 



POWER OF ATTORNEY 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 
Application of; Stuart A. KAUFFMAN et al. 



Application No.; 09/868.981 
Filed: June 22, 2001 

For; AS YSTEM AND METHOD FOR 
THE ANALYSIS AND 
PREDICTION OF ECONOMIC 
MARKETS 



Group Art Unit: To Be Assigned 
Examiner; To Be Assigned 
Attorney Docket No.: 9392-0 16-999 



POWER OF ATTORNEY BY ASSIGNEE 
AED-EXCLUSION OF PfVKNTft K (ft) I™T"F? M C.F.R. 3.71 



Assistant Commissioner for Patents 
Washington, D.C. 20231 

Sir: 

The undersigned assignee of die entire interest in tho above-identified subject 
application hereby appoints: Berj A. Terzian (Reg. No^p060JJDavid Weild, III (Reg. No. 
^21094}Jfpnathan A. Marshall (Reg. N o. 24614), B arry D. Rein (Reg. N o, 22411) , Stanton T. 
Lawrence, III (Reg. N o. 2S736), Ch arles E. McKenney (Reg, No . 22795), P hilip T. Shannon 
(Reg. N o. 24278), F rancis E. Morris (Reg. N o. 246 IS), Ch arles E. Miller (Reg. No. 24576), 
Gidon D. Stern (Reg. No . 27469), John J. Lautar, Jr. (Reg. No . 27814), Bria n M. Poissant 
(Reg. No. 28462) 1 Brian D. Coggio (Reg. No L 27624), Ro ry J. Radding (Reg. No. 28749). 
Stephen J. Harbulak (Reg. N o. 29166), Donald J . Goodcll (Reg. No. 19766), Jame s N. Palik 
(Reg. N o. 2SS10), Th omas E. Friebel (Reg. No . 29258), U ura A, Coruzzi (Reg. No. 30742), 
Jennifer Gordon (Reg. N o. 30753), Gar aldlne F. Baldwin (Reg. No. 31232), Victor nT 
BaWia (Reg. Ko.J^XSw^ B . AbranL3 ^ z N o. 3060S), Stcvcn I. Wallach (Reg. 
No. 35 402), Marc ia H. Sundeen (Rog. No . 30893 ). Paul J. Zeggcr (Reg. N o. 33821), Ed mond 
R. Bannon (Reg. No. 32110), Bruce J. Barker (Reg. No. 33291), Adr iane M. Antler (Reg. No. 
c 32605), T homas G. Rowan (Reg. N o. 34419), Jame s O. Markey (Reg N o. 31636). T homas 
D. Kohlcr (Reg. No327?^^cotip. Stimpson (Reg. No^O^-Gary S. Williams (Reg. 



l 



dci oooesj.l 



POWER OF ATTORNEY 




Famllani (Reg. N o. 35514), Kelly D . Talcolt (Reg. No. 39582), Francis D. Cerrito (Reg. 
No.38100), Anthon y M. Insogna (Reg. No, 35203), Bria n M. Rothery (Reg. No . . 35340), 
Brian D. SifF (Reg. No. 35679). Alan Tencnbaum (Reg. No . 34939), Mic hael J. Lyons (Reg. 
No. 37^386^ Poland T. Stephens (Reg. No. 37,242) and W illiam J. Sipio (Reg, No. 
34 ,5 14), all of Pennie & Edmonds L LP, whose addresses arc U 55 Avenue of the Americas, 
New York, Now York 10Q36, 16 67 K Street N.W., Washington, DC 20006 and 3300 
Hillview Avenue, Palo Alto, CA 94304, all of Pennie Sc Edmonds LLP (PTO Customer No, 
20582), as its attorneys io prosecute this application and to transact all business in the United 
States Patent and Trademark Office connected therewith, said appointment to be to the 
exclusion of the inventors and their attorney^*) in accordance with the provisions of 37 CF.R. 
3.7 1, provided that, if any one of these attorneys ceases being affiliated with the law firm of 
Pennie & Edmonds LLP as partner, counsel, or employee, then the appointment of that 
attorney and all powers derived therefrom shall terminate on the date such attorney ceases 
being so affiliated. 

An assignment of the entire interest In the above-identified subject application: 

[ ] was recorded on at reel/frame / . 

[X] is submitted herewith for recording. 

Please dirtict all correspondence for this application to customer no. 20582, 



ASSIGNEE: 



Signature: 
Typed Name: 




Position/Title: 



Robert R. MacPonald 
IfVesidcnt 



Address; 



317 Pasco dc Peralta 



Santa Fc. Ne w Mexico 87501 



Dale: 
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DECLARATION FOR NON-PROVISIONAL PATENT APPLICATION* 



As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below at 201 et seq. beneath my name. 

I believe I am the original, first and sole inventor if only one name is listed at 201 below, or an original, first and joint inventor if plural names 
are listed at 201 et seq. below, of the subject matter which is claimed and for which a patent is sought on the invention entitled 

A SYSTEM AND METHOD FOR THE ANALYSIS AND PREDICTION OF ECONOMIC MARKETS 

and for which a patent application: 

□ is attached hereto and includes amendment(s) filed on ofappucabu) 

h was filed in the United States on June 22, 2001 as Application No. 09/868,981 (for *ci ara non no< acco^^ng appil cano nJ 
with amendment(s) filed on o/ applicable) 

a was filed as PCT international Application No. PCT/US99/30640 on December 22, 1999 and was amended under PCT Article 19 on M 

applicable) ' J 

I hereby state that I have reviewed and understand the contents of the above identified application, including the claims as amended by anv 
amendment referred to above. * 

I acknowledge the duty to disclose information known to me to be material to patentability as defined in Title 37, Code of Federal Regulations 
§1.56. 

I hereby claim foreign priority benefits-nnder Title 35, United States Code, §1 19(a)-(d) of any foreign application(s) for patent or inventor's 
certificate listed below and have also identified below any foreign application for patent or inventor's certificate having a filing date before that 
of the application on which priority is claimed: 



EARLIEST FOREIGN APPLICATION(S), IF ANY, FILED PRIOR TO THE FILING DATE OF THE APPLICATION 


APPLICATION NUMBER 


COUNTRY 


DATE OF FILING 
(day, month, year) 


PRIORITY 
CLAIMED 








YES □ NO □ 








YES □ NO □ 








YES □ NO a 



I hereby claim the benefit under Title 35, United States Code, § 1 19(e) of any United States provisional application(s) listed below. 



PROVISIONAL APPLICATION NUMBER 


FILING DATE 


60/113,502 


December 22, 1998 











I hereby claim the benefit under Title 35, United States Code, §120 of any United States application(s) listed below and insofar as the subject 
matter of each of the claims of this application is not disclosed in the prior United States application in the manner provided by the first 
paragraph of Title 35, United States Code §1 12, 1 acknowledge the duty to disclose information known to me which is material to patentability 
as defined in Title 37, Code of Federal Regulations, §1.56 which became available between the filing date of the prior application and the 
national or PCT international filing date of this application: 



NON-PROVISIONAL 
APPLICATION SERIAL NO. 


FILING DATE 


STATUS 


PATENTED 


PENDING 


ABANDONED 


PCT/US99/30640 


December 22, 1999 




X 

























* for use only when the application is assigned to a company, partnership or other organization. 

(I) ' 



DCI -300946.1 



I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are 
believed to be true; and further that these statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements 
may jeopardize the validity of the application or any patent issuing thereon. 


2 
0 
1 


FULL NAME 
OF INVENTOR 


LAST NAME 

KAUFFMAN 


FIRST NAME 

Stuart 


MIDDLE NAME 

A. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa_Fe — 


STATE OR FOREIGN COUNTRY 

NM N<V\ 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

1 8 1 1 S. Caimno Cruz Blanca 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE C/FJtWENTOR 20 1 ^t^^^^^^^^ 


DATE 


2 
0 
2 


FULL NAME 
OF INVENTOR , 


/ 

LAST NAME 

MAC READY 


FIRST NAME 

William 


MIDDLE NAME 

G. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 

Canada 


POST OFFICE 
ADDRESS 


""STREET 1 ' 

339!4Delgado 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 202 


DATE 


2 
0 
3 


FULL NAME 
OF INVENTOR 


LAST NAME 

PALMER 


FIRST NAME 

Richard 


MIDDLE NAME 

G. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Durham 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 

United Kingdom 


POST OFFICE 
ADDRESS 


STREET 

4030 Forestdale Drive 


CITY 

Durham 


STATE OR COUNTRY 

NC 


ZIP CODE 

27712 






SIGNATURE OF INVENTOR 203 


DATE 

10-2 - -ZOO J 


2 
0 
4 


r ULL NAJvlb 
OF INVENTOR 


LAST NAME 

COWAN 


FIRST NAME 

Jack 


MIDDLE NAME 

D. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Chicago 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

5448 Ridgewood Court 


CITY 

Chicago 


STATE OR COUNTRY 
IL 


ZIP CODE 

60615 






SIGNATURE OF INVENTOR 204 


DATE 


2 
0 
5 


FULL NAME 
OF INVENTOR 


LAST NAME 

TEILTELBAUM 


FIRST NAME 

Daniel 

s 


MIDDLE NAME 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 

NM 


COUNTRY,OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

618 Old Santa Fe Trail #G 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 205 


DATE 



(2) 



DO - 300946.1 



( hereby declare that alt statements made herein of my own knowledge arc true and tfcal all sUtttncnti made oh Jn/orrnation and belief arc 
believed to be iTM«i and tiinhfir iKit th«c siatcmepis were raid* with the knowledge thai wlllflil <aUc stataoientt and the Hke so mode vc 
punlaihabla by fine or imprisonment or both, under Section )00l oifThlc 18 of the Unltad State) Code and thai such willftf f&Uc sutamenU 
may jeopardize the validity of the application or any patent issuing thereon. 





FULL NAME 
OP INVENTOR 


Last wami 
OUTXIN 


Aicxandor 




2 
0 
6 


RESIDENCE & 
■ CITIZENSHIP 


crrv 

Santa Fc 


TTATt Oft COUNTRY /A A 


eotfttw of cmnuwtir 




post omc£ 

ADDRESS 


rtJiticr 

\$48 AvI<Ja dc Las Americas 


CITY 

Santa Fe 


mte oixouwntY 
NM 


a* coos ■ 
37501 






IfOUATVftE <U iKYfNTO* JOI 






FULL NAME 
Of INVENTOR 




Bruce 


HIDOUWAMf 

K 


7 


! : -Mi353ENCE^ ; '^ 

crnzfcNSHH* 


ijcttv ■. . 
_Sftij3.Ee— 


•waWo* kjajw cowraV . . ' 


' United Stain 




K>ST OFFICE 
ADDRESS 


muurr 
102$ Don Cubaro 


crrv 1 
Santa Fe 


jtati o« camrrty 
NM 


87501 








DATE 

t/ 24/ £00 / 




FULL NAME 
OP INVENTOR. 


MAKHANKQV 


Vladimir 


Q. 


2, 
0 
8 


RESIDENCES 
CITIZENSHIP 


crfY 




COUNTRY Of CTITOSN5 Wf 

Russia 




POST OFFICE 
ADDRESS 


3TJUUT 

4143 Chapatton Plaoi 


crrr 
Santa Fe 


tT<rsoncouwTi.v 
NM 


Xlf COM! 

87301 








-) ^ 


&ATS 




FULL NAME 
OF INVENTOR 


DARLEY 


mjutt w*vc 
^Vtiicc 


HTOOit NAME 

M. 


2 

0 


R£Sp£NC&£ _ 


crrv 


HATE Ot POU&QN 


covfrriLY o p aru»r*swp 






TOSTOFFJCE 
ADDR£$S 


794 Caile Altimira 


crrv 

SaauF* 


mTIo*co^W^*y , 
NM 


to>cooe 
8730) 










DATE ^ 






FULL NAME 
OF INVENTOR 








2 
0 
10 


RESIDENCE <fe 
CITIZENSHIP 




IT A« OIL COUNT*. T 


tfOLVT* Y Of OTtEIfclHl* 




POST OFFICE 
ADDRESS 




crrv 


ETAT* OR COvVTRV 






WW»TUM Of fNYXNTOU 33i 


D»Tt 





T htireby declare iha* all statements nude herein of o\y own knowledge axe true and thai al] tuwmcatj rtwk oa iaformadon and belief urc 
believed to be vrue, and further lh»t these jc^cmcftts wera raid* w\\h tfic knowtodflc thai Willfbl f«J jc ff*f/tnv5n« And Th<3 (ike so mode uc 
pun^tublc by fine or ImpnwamOTiL or bovh. imdtjr S4ctfon 1001 of Title 1 8 of the United Stated Code tad th«t such willM falic tuumenu 
may jcop*rdi2c the validity of the application or any paLtut issuing thereon. 


2 
0 
6 


FULL NAME 
OFIKVEKTOR 


OUTXJN 






RESIDENCE & 
■ CITIZENSHIP 


cmr 

Santa Fc 


rrjsTw c* ;o*a on country 
NM 


cotftftHY of crrnawjMir 


POST OFFICE 
ADDRESS 


1 545 Avida dc Las Americas 


Santa Fe 


mrta*coufTrv 
NM 


37501 








OA. TIC 


7 


FULL NAME 
Of INVENTOR 


SAWHILL 


Bruce 


MlDOUKAMf 




jarv . 
Santa Fe 


^rAH 0*. fO*i*3W COUKTlilf 


'* United Stall* 


tOST OFFICE 
ADDRESS 


1023 Don Cutwro 


cmr ' 
StntaFe 


itati o* cocvrtY j iff coot 
NM j R7501 






ticKATtAi of MVfxva* m , — . 


HATS. 


0 
6 


FULL NAME 
OP INVENTOR 


MAKKANKOV 


Vladimir 


kingu: hAwt 
G. 


RESIDENCE <1 
CITIZENSHIP 


crrv 

Santa Fe 


IT ATI OR fO»£JGN COX/NT* V 

NM 


coUKYHYOf CmasKSWf 
Russia 


POST OFFICE 
AUORH5S 


4 J 48 Chaperron Plico 


cm 

Santa Fe 


nATCOnCOWTC*' 

NM 


ii/ cone 
$7301 




s (oua rvtcxx^^^ • 


fUTX 


2 
0 


FULL NAME 
OF INVENTOR 


DARLEY 


Vlncc 


M. 




crrr 

TT ; — "i — i > ■ '•■ 


tTATf 0 1 tCXJiQH C(kJT<T».Y 


cog^rm or aTU*<s*iF 




TOST OFFICE 
ADDRfeSS 


794 Callc AJtunirt 


■crrv 


irA^O^CCfWfKY* j fc»COO< 

NM J t7i0J 








A * 

CWLTt 

(( jto J-z^) 




FULL NAME 
OF INVENTOR 








2 
0 
10 


RESIDENCE & 
CmZENSHTP 


crrv 


rune on. poiuboncquktut 


< ocvT* r Of cnxpD/i fi 1 ^ 


POST OFFICE 
ADDRESS 


'iTitm 


OTY 


trxTt o*.cout/T>v 


n* coot 




WOW*.T\J»UJ Of mtNToft 33i 


t)4tt 





1 licrcby declare lhal ail statements made herein of my own knowledge are true and that a! I statements made on information aiid belief are 
believed to be true; and further that these statements were made with the knowledge that willful false statements and the like so made arc 
punishable by Fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements 
may jeopardize the validity of the application or any patent issuing thereon. 



2 
0 
6 


FULL NAME 
OP INVENTOR 


LAST NAME 

OUTKIN 


FIRST NAME! 

Alexander 


MIDDLE NAME 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fc 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY Of CITIZENSHIP 

Russia 


POST OFFICE 
ADDRESS 


STREET 

1548 Avida de Las Americas 


CITY 

Santa Fc 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 201 


DATE 


2 
0 
7 


FULL NAME 
OF INVENTOR 


LAST NAME 

SAWHILL 


FIRST NAME 

Bruce 


MIDDLE NAME 

K. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 

NM 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


sTRcer- 
1028 Don Cubero 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 




SIGNATURE OF INVENTOR 202 y n. 


DATE 


2 
0 
8 


FULL NAME 
OF INVENTOR 


LAST NAME 

MAKHANKOV 


FIRST NAME 

Vladimir 


MIDDLE NAME 

G. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

Russia 


POST OFFICE 
ADDRESS 


STREET 

4148 Chaparron Place 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 




SIGNATURE OF INVENTOR 203 


DATE 




2 
0 
9 


FULL NAME 
OF INVENTOR 


LAST NAME 

DARLEY 


FIRST NAME 

Vince 


MIDDLE NAME 
M. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

United Kingdom 


POST OFFICE 
ADDRESS 


STREET 

794 Callc Altamira 


CITY 

Santa Fe 


STATE OR COUNTRY 

Nrvt 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 2<M {^^^^^^ 


DATE 

(t JiO ) 


2 
0 
10 


FULL NAME 
OF INVENTOR 


LAST NAME 


first Name 


MIDDLE NAME 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 


POST OFFICE 
ADDRESS 


STREET 


CITY 


STATE OR COUNTRY 


ZIP CODE 






SIGNATURE Of INVENTOR 205 


DATE 



I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are 
believed to be true;. and furtner that these statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statements 
may jeopardize the validity of the application or any patent issuing thereon. 


2 
0 
I 


FULL NAME 
OF INVENTOR 


LAST NAME 

KAUFFMAN 


FIRST NAME 

Stuart 


MIDDLE NAME 
A. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

181 1 S. CaFrrjno Cruz Blanca 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 






^signat^^ 


DATE 


2 
0 
2 


FULL NAME 
OF INVENTOR 


s 

LAST NAME 

MACREADY 


FIRST NAME 

William 


MIDDLE NAME 

G. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

Canada 


POST OFFICE 
ADDRESS 


STREET 

339 l /r*>elgado 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 202 


DATE 


2 
0 
3 


FT II I NAMF 
OF INVENTOR 


LAST NAME 

PALMER 


FIRST NAME 

Richard 


MIDDLE NAME 

G. 


RESIDENCF & 

IV! — t i^J X XS L/ i T *w X\j 

CITIZENSHIP 


CITY 

Durham 


STATE OR FOREIGN COUNTRY 

NC 


COUNTRY OF CITIZENSHIP 

United Kingdom 


POST OFFICE 
ADDRESS 


STREET 

4030 Forestdale Drive 


CITY 

Durham 


STATE OR COUNTRY 

NC 


ZIP CODE 

27712 






SIGNATURE OF INVENTOR 203 


DATE 


2 
0 
4 


FULL NAME 
OF INVENTOR 


LAST NAME 

COWAN 


FIRST NAME 

Jack 


MIDDLE NAME 

D. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Chicago 


STATE OR FOREIGN COUNTRY 

IL 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

5448 Ridgewood Court 


CITY 

Chicago 


STATE OR COUNTRY 

IL 


ZIP CODE 

60615 




SIGNATURE OF INVENTOR 204 


T 

□ATE X 


2 
0 
5 


FULL NAME 
OF INVENTOR 


LASTNAME fT^U^t^t^ 

■ TE^TELDAUM- > ^ 


FIRST NAME 

Daniel 


MIDDLE NAME 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 

Unitetl States 




POST OFFICE 
ADDRESS 


street V©6/ AJ.^ i±f^|<-| 


CITY ftfijLit*tT<^ 


STATE OR COUNTRY 


ZIP CODE 


<6 1 0 Old Gain* Ti TiaTl #G» 


•87501-23^0^ 




SIGNATURE OF INVENTOR 20> 

/ij^-^f /c\job ■ 


DATE 

(o • Of • 2^> 1 



DCl - 300946.1 



I hereby declare that all statements made herein of my own knowledge arc true and that all statements made on information and belief are 
believed to be true; and further that these statements were made with the knowledge that willful false statements and the like so made arc 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United Slates Code and that such willful false statements 
may jeopardize the validity of the application or any patent issuing thereon. 





FULL NAME 
OF INVENTOR 


LAST NAME 

KAUFFMAN 


FIRST NAME 

Stuart 


MIDDLE NAME 

A. 


2 
0 
1 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR. FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

United States 




POST OFFICE 
ADDRESS 


STREET 

1811 S. Ca^Frrjno Cruz Blanca 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATTJRE C/FflWENTOR 201 


2^ 


DATE 






-r* 










FULL NAME 
OF INVENTOR 


LAST NAME 

MAC READY 


FIRST NAME 

William 


MIDDLE NAME 

G. 


2 
0 
2 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 

NM 


COUNTRY OF CITIZENSHIP 

Canada 




POST OFFICE 
ADDRESS 


STREET 

339Vf5elgado 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 






SIGNATURE OF INVENTOR 202 


DATE 






FULL NAME 
OF INVENTOR 


LAST NAME 

PALMER 


FIRST NAME 

Richard 


MIDDLE NAME 


2 
0 
3 


RESIDENCE & 
CITIZENSHIP 


CITY 

Durham 


STATE OR FOREIGN COUNTRY 
NC 


COUNTRY OF CITIZENSHIP 

United Kingdom 




POST OFFICE 
ADDRESS 


STREET 

4030 Forestdale Drive 


CITY 

Durham 


STATE OR COUNTRY 

NC 


ZIP CODE 

27712 






SIGNATURE OF INVENTOR 203 




DATE 












10-2. - Zoo 1 




FULL NAME 
OF INVENTOR 


LAST NAME 
COWAN 


FIRST NAME 

Jack 


MIDDLE NAME 

D. 


2 
0 
4 


RESIDENCE & 
CITIZENSHIP 


crry 
Chicago 


STATE OR FOREIGN COUNTRY 

IL 


COUNTRY OF CITIZENSHIP 

United States 




POST OFFICE 
ADDRESS 


STREET 

5448 Ridgewood Court 


CITY 

Chicago 


STATE OR COUNTRY 
IL 


ZIP CODE 

60615 






SIGNATURE OF INVENTOR 204 


* * 

DATE 






FULL NAME 
OF rNVENTOR 


LAST NAME 

TEILTELBAUM 


FIRST NAME 

Daniel 


MIDDLE NAME 


2 
0 
5 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY.OF CITIZENSHIP 

United States 




POST OFFICE 
ADDRESS 


STREET 

618 Old Santa Fe Trail #G 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 




SIGNATURE OF INVENTOR 205 


DATE 





IX' J - .10U'M6 1 



[hereby declare that *U statement* made herein of my own knowledge trc inic and thai all lUtcmcnti mtdc on infonnilioo md behef arc 
believed to be mjc; and ftffrficr ihai th«e statements wtrc midc with the knowledge ihti willful false it*icmcnta and ihc like so made an. 
punl^ablc-byflncorlrapriwnmcni, or boiti, under Section 1001 ofTi:le lg of the United Sta las Code and iha; cuch wlllftjl (alac ctttemencs 
m&yJ*opardlt« the validity af the application or &ny pal«nt issuing (her ton. 



2 
0 
1 


FULL NAME 
OF INVENTOR 


kauffman 


Stu*n 


WtDOLCMAWC 

A. 


RESIDENCE A 
CITIZENSHIP 


C(TY 

Sdnm Fc 


ITXTT0* COUkTHY 

NM 


COVNTKYOf CTTTZZ>nwr 

Untied Suies 


POST OFFICE 
ADpHES$ 


1811 S. Curtw Cruz Blanca 


cttt 
y Santa. Fe 


STATE OC COUXTKY 

NM 


w cooe 
87501 








DAT* 


2 
0 
2 


■ ■ r-< 

FULL NAME 
OF INVENTOR 


Last namm 

macready 


William 


a 


RESIDENCE & 
CITIZENSHIP 


CrtY 

S*nta Fa 


CTATS C* POIIEICM eot*rpLY 

NM 


COtfKTfLY Of CTTVUXtHlr 

Canada 


POST OFFICE 
ADDRESS 




Santa Fe 


itati o*. datJwmY 
NM 


BP C0OC 

rotrr 






hqMatuk£ of iwvorroit. *o 




2 
0 
3 


FULL NAME 
OF INVENTOR 


LASTNAXl 

PALMER 


Richard 


MlDgUTHAMS: 

G. 


RESIDENCE Sc 
CTITZENSHIP 


CITY 

Durham 


ctati ot roaricw coomnr 
NC 


covkiily or araiww 
United Kingdom 


POST OFFICE 
ADDRESS 


xruurr 

4030 Forcytdilc Drive 


crrr 

Durham 


state o«. count* y 
NC 


Or coot 
27712 






*K3<atv«£ 0» rwvarroa i«i 


OA Til 


2 
0 
4 


FULL NAME 
OF INVENTOR 


L*JT «AmS 

cowan 


CUTTKAXC 

Jack 


KiooccvAxa 
D. 


RESIDENCE «. 
.CITIZENSHIP 


CITY 

Chicago 


iTAte o*i rcwziow country 


COWTHY 0/ CTTCE^tUXP 

United States 


POST OFFICE 
ADDRESS . 


STTUVT 

5448 RJdgcwood Court 


cmr 
Chicago 


JTAT«ORCOfUWnt.Y 

IL 


2ft COOG 

60615 






SlC*mUR£ or IKVEMTOl 2<M 


DATE 

* * 


2 
0 
5 


FULL NAME 
OF INVENTOR 


^TEyrciAAUM ^ 


Daniel 


mooCCUAMK 


RESIDENCE <fe 
CITIZENSHIP 


cmr 


rr*Tr oc »o«uuow oournrT 


cxKJKmYor ccnzswHir 
United Sttiu 


POST OFFICE 
ADDRESS 


^10 0( d CmiLu Pu'Tl Cfl flG* 


Fo 


STaTE O*. C OCA/T* V 


rrcoac 






M^r^f MflA 


a* re 



I hereby declare that all-statements made herein of my own knowledge arc true and thai all statements made on information and belief 
believed to be true; and further that these statements were made with the knowledge that willful false statements and the like so made i 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such willful false statei 
may jeopardize the validity of the application or any patent issuing thereon. 



2 
0 
I 


FULL NAME 
OF INVENTOR 


LAST NAME 

KAUFEMAN 


FIRST NAME 

Stuart 


MIDDLE NAME 

A. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATF OR FOREIGN COUNTRY 

NM 


COUNTRY OF CITIZENSHIP 

United States 


POST OFFICE 
ADDRESS 


STREET 

1811 S. O^rrrjno Cruz Blanca 


CITY 

Santa Fe 


STATE OR COUNTRY 
NM 


ZIP CODE 

87501 




SIGNATURE C/f^VENTOR 201 ^^^^^J^^ 


DATE 

9 'Z 




2 
0 
2 


FULL NAME 
OF rNVENTOR 


LAST NAME 

MACREADY 


FIRST NAME 

William 


MIDDLE NAME 

G. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Santa Fe 


STATE OR FOREIGN COUNTRY 
NM 


COUNTRY OF CITIZENSHIP 

Canada 


POST OFFICE 
ADDRESS 


STREET 

339 ! /z-©elgado 


CITY 

Santa Fe 


STATE OR COUNTRY 

NM 


ZIP CODE 

87501 




SIGNATURE OF INVENTOR 202 


DATE 




2 
0 
3 


FULL NAME 
OF INVENTOR 


LAST NAME 

PALMER 


FIRST NAME 

Richard 


MIDDLE NAME 

G. 


RESIDENCE & 
CITIZENSHIP 


CITY 

Durham 


STATE OR FOREIGN COUNTRY 

NC 


COUNTRY OF CITIZENSHIP 

United Kingdom 


POST OFFICE 
ADDRESS 


STREET 

4030 Forestdale Drive 


CITY 

Durham 


STATE OR COUNTRY 

NC 


ZIP CODE 

27712 




SIGNATURE OF INVENTOR 203 


DATE 




2 

V 

4 


pr ir r XT A K/tV 

OF INVENTOR 


LAST NAME 

COWAN 


FIRST NAME 

Jack 


MIDDLE NAME 

D. 


CITIZENSHIP 


CITY 

Chicago 


STATE OR FOREIGN COUNTRY 

IL 


COUNTRY OF CITIZENSHIP 

United States 


ir VJo 1 UrrlCb 
ADDRESS 


STREET 

5448 Ridgewood Court 


CITY 

Chicago 


STATE OR COUNTRY 

IL 


ZIP CODE 

60615 






SIGNATURE OF INVENTOR 204 


» 1 

DATE * 






FULL NAME 
OF INVENTOR 


LAST NAME |T6t/^Wt^\ 

• XElffTELBAUM- p ^ 


FIRST NAME 

Daniel 


MIDDLE NAME 


2 
0 

5 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 

-NM \ff\ 


COUNTRY OF CITIZENSHIP 

Unitetj States 




POST OFFICE 
ADDRESS 


STREET H00i AMtt # f^|<-{ 


CTY Pifk^^CV^ 


STATE OR COUNTRY 


ZIP CODE 




SIGNATURE OF INVENTOR 20> 


DATE 

(o • 0<f • 2A° 1 



